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Chapter 1 


Introduction 

1.1 Background and Motivation 

The flight safety of an airplane critically relies upon the proper operation of various 
communication, weather, navigation and other avionic systems [1]. Many of these 
systems have antennas and/or other structures on the skin of the aircraft. The 
operational characteristics of any antenna on the airplane include the effect of the 
various structures making up the airframe and the electromagnetic coupling among 
the multitude of antennas and transponders. Since antennas are usually designed 
assuming that they are isolated and mounted on flat ground planes, the actual 
pattern/performance of the antenna can significantly differ from their theoretical 
designs. Hence, for the safe operation of the aircraft, it is vital to predict the 
electromagnetic performance of each system in the presence of others. 

With the increasing number of airplanes in the sky, avoidance of mid-air collisions 
has assumed new importance [2]. To this end, the Federal Aviation Administration 
has mandated that all aircraft that can carry more than 30 passengers be fitted 
with an advanced Traffic-alert Collision Avoidance System (TCAS) [3]. The TCAS 
is basically a direction finding and tracking radar, that determines the bearing of 
a nearby aircraft from a signal received from that aircraft. There are a number of 
tracking radars described in the literature; however, due to cost, space, weight as 
well as other constraints, only a few of these can be used in airborne applications [4]. 
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The basic operation of the TCAS is as follows: The TCAS equipped airplane 
transmits a directional query signal. This query triggers a transponder on nearby air- 
planes, if so equipped. The digitally coded, transponder return signal is received by 
the TCAS and processed to detect the presence of an intruder. It may also estimate 
the bearing, velocity etc., of the other airplane. For convenience, the TCAS equipped 
aircraft will be referred to as the own airplane and the transponder equipped nearby 
airplane will be called the intruder . 

Presently, there are three generations of the collision avoidance systems, named 
TCAS I, II and III respectively. The TCAS I system, the oldest, could only inform 
the pilot about the presence of other aircraft in his vicinity. It could not advise 
the pilot of evasive maneuvers to take to avoid a potential collision. The TCAS II 
system, currently in the production phase, has conflict resolution logic, which can 
suggest that the pilot should climb or dive [5]. The most recent system, the TCAS 
III has the ability to determine the bearing of the intruder aircraft accurately, and 
can issue to the pilot an advisory to bank left or right in addition to climb or dive. 
This is possible because, along with the bearing, the TCAS III system can also 
determine the altitude, velocity, and distance of the intruder from the transponder 
response. More information about these systems can be found in [6], [7], [8]. 

There are many factors that adversely affect the performance of the TCAS in 
estimating the bearings of surrounding aircraft. Even though it is not very hard 
to predict the bearing of an intruder under ideal design conditions, it is much more 
difficult to predict the performance of the same system when the antenna is mounted 
on an actual aircraft. This is due to the distortion of the radiation pattern by the 
curvature of the fuselage, scattering by the wings, engines, horizontal and vertical 
stabilizers as well as other nearby antennas located on the skin. The cost of making 
full scale measurements and evaluating these systems under actual flight conditions is 
very high, time consuming, and error prone. Therefore, there is a need for alternative 
methods to study the performance of various systems. 
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One common method is to nse scale models to make measurements [9]. This is 
inflexible in that one can not study all possible locations for a TCAS antenna on 
the fuselage of various aircraft. Another possibility is to use computer models to 
calculate the radiation patterns in the near or far field, and process the signal on 
a computer. The latter method is adopted here, and in the succeeding pages, the 
TCAS system performance is evaluated under various adverse conditions. 

1.2 Synopsis 

In this report, two different TCAS antennas are studied in detail. They are classified 
based on the operation. These are: 

1. Comparison of Relative Amplitude system (CRA). 

2. Spiral Phase Antenna system (SPA), and 

The antennas for both systems consist of a four element circular array of monopoles; 
however, their operation and their performance in the presence of nearby scatterers 
is quite different. The CRA TCAS finds the bearing by comparing the amplitudes 
of the received signals; whereas, the SPA TCAS determines the bearing by using 
only the phase information contained in them. In this work, three main topics are 
covered: 

• Modeling isolated TCAS antennas on aircraft. 

• Simulation of TCAS performance in the presence of nearby antennas. 

• Modeling of engine to better account for its effect on the TCAS. 

Chapter 2 introduces the two systems. Their operation, modeling the antennas 
and the issues involved in modeling their performance are discussed. Models for the 
CRA and the SPA TCAS systems are developed. These are used in conjunction 
with the high frequency models of various airplanes. The radiation patterns are 
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obtained using the airborne antenna code, called NEWAIR3, based on the uniform 
geometric theory of diffraction [10], [11], [12]. The impact of various scatterers on the 
TCAS is characterized by the bearing error, defined as the difference between the 
system computed bearing and the actual bearing of an intruder. Simulation of the 
TCAS and computation of the bearing over all azimuth results in a bearing error 
curve. Under ideal conditions, the bearing error is zero. Due to the scatterers and 
other antennas on the airframe, which distort the received signals, this is not the 
case. The magnitude of the error indicates the relative effect of various scatterers 
on the TCAS [13], [14]. The salient features and properties of the error curves are 
also discussed in Chapter 2. 

In Chapter 3, models are derived for different airframes and error curves are 
computed for each TCAS antenna at various positions on the fuselage. The error 
curves are then used to compute the standard deviation of the bearing error in a 
specified angular sector, which is used as a figure of merit to evaluate the performance 
of the TCAS at that location. It can be either used to compare the two systems 
TCAS on the same fuselage at two different locations or the same system on different 
airframes at similar locations. Three airplanes were chosen: Boeing models 737-300, 
747-200 and 727-200. The variation of the standard deviation of the bearing error as 
a function of distance from the nose are presented. The curves of standard deviation 
versus the location on the fuselage are provided for both the top and bottom mounted 
SPA and CRA antenna systems [15]. It is shown that the various aircraft scatterers 
play an important role in the determination of the angle of arrival (AO A) errors. 
The need for a better engine model is demonstrated for some cases, specifically when 
an engine is directly illuminated by the TCAS. 

The effect of two nearby antennas mounted aft of the TCAS antenna at optimum 
location on the behavior of the TCAS is considered next in Chapter 4. The scattering 
from the blade antennas in the vicinity were included in the patterns of the TCAS 
antenna via superposition [16]. Two representative antennas were chosen for this 
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study; one a VHF blade, and the other an L-band blade. Modeling methods for 
these two antennas are different and hence are discussed separately. First, using a 
moment method analysis, the scattering fields were computed [17]. Next, using an 
antenna synthesis procedure based on the method of least squares, models for the 
blade antennas were developed. It is found that these nearby antennas have a bigger 
impact on the SPA TCAS than the CRA system. This is especially true in the 
forward directions. The reasons for this phenomenon are discussed, citing examples 

from the study on the Boeing 737-300 [15]. 

Chapter 5 discusses the modeling of the aircraft engine as a circular waveguide 
cavity. The waveguide was terminated using a perfectly conducting flat plate. The 
patterns of the scattering from the engine were obtained via modal techniques and 
then superimposed on the patterns of the TCAS [18], [19], [20]. It was discovered that 
the engine model consisting of flat plates overestimates the error in certain regions 
and under-estimates in others. The waveguide model leads to a more uniform bearing 
error. It is also seen that the CRA is affected much less than the SPA system. The 
performance of the system can be improved slightly by moving the antenna forward. 
The advantages and disadvantages of modeling the engine as a terminated waveguide 

versus a closed rectangular box are discussed. 

The results are summarized in Chapter 6. Guidelines are devised to locate the 
TCAS antenna on the fuselage without compromising the performance in the pres- 
ence of various scattered. The results indicate that for best performance, one should 
place the TCAS antenna towards the forward portion of the fuselage. For forward 
look angles, it is found that the CRA TCAS performs better than the SPA system 
under all situations considered in this report. Consequently, installation of the CRA 
antenna poses fewer difficulties than the SPA antenna. This could be an important 
factor when the space available on the fuselage is small or there are other nearby 
scattered. Based on results obtained so far, antenna placement on other aircraft is 
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briefly considered. Finally, some techniques that, might improve the performance of 
the TCAS antennas in general are suggested. 
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Chapter 2 

Operation of TCAS 


2.1 Introduction 

Two angle of arrival systems are introduced in this chapter. They are classified as: 

• Comparison of Relative Amplitude (CRA) system. The CRA TCAS operates 
by comparing the amplitudes of the received signals via different beams. The 
phase information is normally ignored. 

• Spiral Phase Antenna (SPA) system. A competing TCAS system, that finds 
the bearing of the intruder from the phase of the received signals. This algo- 
rithm usually discards the amplitude information. 

The TCAS antennas are modeled and the systems simulated on a digital computer. 
The properties of the the two systems are explored using the bearing error transfer 
functions. Some bearing error curves for the CRA TCAS are found to have disconti- 
nuities at some specific angles. It is discovered that the SPA TCAS is more sensitive 
to nearby scatterers than the CRA system. 

2.2 Comparison of Relative Amplitude System 

The CRA TCAS employs a four element circular array. This array transmits a beam 
in one of four selectable directions. The array is mounted on the fuselage and aligned 
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Figure 2.1: Typical directional interrogation beam for the amplitude TCAS. 

in such a way that the beams are directed in the 0°, 90°, 180° and 270° of relative 
bearing with respect to the forward direction. The beams have a 3-dB beamwidth 
of 90° ± 10° for all elevation angles between +20° and —15°. Note that 0° elevation 
corresponds to the standard azimuth pattern cut. A sketch of the desired receiver 
pattern is shown in Figure 2.1. 

The frequency of transmission of the query signal is 1030 MHz. The received 
signal is received individually by each element of this antenna array at 1090 MHz. 
The magnitudes of the four received signals are then compared and the difference 
between the largest two computed. Comparison of this difference signal with similar 
values on a calibration curve yields the angle of arrival (estimated bearing of the 
intruder). 

The TCAS antenna is modeled by four short z-directed monopoles on a flat 
ground plane as depicted in Figure 2.2. The current distributions on the monopoles 
are assumed to be sinusoidal. Let the excitation of the individual elements be given 
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CRA antenna 

Figure 2.2: Geometry of the 4 element CRA TCAS model. 


b y 

/, = A,e ja ' (2.1) 

where, I, is the complex excitation, A, is the amplitude, and a, is the phase of element 
i, where i = 0, 1,2,3, and j = y/—l. The elements are numbered counterclockwise 
from the X-axis as shown. 

The far field antenna pattern Eg(0,<f>) is the product of the element pattern 
F{r,6,<f>) and the array factor AF(6,<f>) such that 


E e (0,<f>) = F(r,e,t)-AF(0,t) 


( 2 . 2 ) 


where, the azimuth angle <f> is measured counterclockwise from the forward direction 
(towards the nose), and 6 is measured from the vertical as shown in Figure 2.3. 

The element pattern of a dipole of length I or a monopole of length 1 / 2 on a 
ground plane can be written as [21] 


F( r ,e,<j>) = 


jqe J * r 

2trr 



(2.3) 
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Figure 2.3: The aircraft coordinate system. Azimuth angle is measured counter- 
clockwise from the forward direction. Elevation angle is measured from the horizon. 


and the array factor is given by 

AF{0 y <f>) = E i4 1 -c' (Qi+fca " n ® co,( *"* )) ( 2 - 4 ) 

1=0 

where, a is the radius of the array, and k is the free space wavenumber, and k = 2tt/A. 
^, = if is the angular position of the i lh element, 7/ is the impedance of free space, 
and (r, 6, <f>) are the standard spherical coordinates. 

Since the number of elements is small, the amplitudes are set to a constant value 
A, = A, and only the phase of the excitation is modified so as to get a maximum in 
any one of the four directions given by fa = j\,j = 0,1, 2, 3, corresponding to nose, 
port side (left), tail, and right (starboard), respectively. Also, in order to achieve 
a symmetric azimuthal pattern, the two elements perpendicular to the axis of the 
beam heading must have the same phase. From these conditions, for example, for 
a beam maximum in the direction of the nose (6 = 0 v ,<f> = 0°), one possibility is to 

pick 

<*o = arbitrary (2.5a) 
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at] = a () + ka sin 9 t) 


(2.5b) 

a 2 — a 0 -f 2ka sin 

and 

(2.5c) 

a 3 - ai. 


(2.5d) 


The weights were calculated for each beam heading <f>j y and = 80° , because this 
0u corresponds to the elevation angle of most interest, namely 10°. These weights 
were then kept fixed in further calculations. To scan the space, one simply adjusts 
the phase distribution of the antenna elements. 

The radiation patterns of the TCAS antenna in this synthesis, as well as 
other simulations were obtained by using a FORTRAN program called NEWAIRS. 
NEWAIR3 is based on the high frequency method, the Uniform Geometric Theory 
of Diffraction (UTD). This is sometimes referred to as the aircraft code or the air- 
borne antenna code. The code can be used to compute the radiation patterns of 
short monopoles and slots on a curved body. The body is modeled by a composite 
ellipsoid, joined smoothly at the location of the antenna. Other structural features 
such as wings, tail etc. are simulated by perfectly conducting (PEC) flat plates. 
Details may be found in the literature [12, 11]. 

From a number of simulations with varying parameters of the array, a radius of 
a = 2.28" and antenna length of Z = A/4 were chosen at the frequency of simulation 
1090 MHz. The length of the antenna elements affects only the elevation plane 
patterns, and the effect is small for l < A/4. As an example, the four beams were 
calculated on the fuselage of a Boeing 727-200 at a distance of 380" from the nose. 
The computer model consisted of only the fuselage. The beams are depicted in 
Figure 2.4 and the important properties are tabulated in Table 2.1. It is noted that 
the radiation characteristics are different for all four beams because of the effect of 
fuselage curvature. 

The model for the operation of the amplitude system proceeds by first calculating 
the calibration or lookup table. For the CRA system, this consists of the difference 
of the dB magnitude of the received signals for the antenna mounted on the fuselage 
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Figure 2.4: Four beams for top mounted CRA array located 380" from the nose on 
the fuselage of Boeing 727-200. 


Table 2.1: Parameters of the four CRA beams 380" from the nose on Boeing 727 
fuselage (radius = 2.28"). 


Beam 

Bearing 

Beam width (°) 

Drop at 
90° (dB) 

3 dB 

10 dB 

Ideal 

00 

180 

10.00 

Nose 

85.1 

160.5 

13.2 

Port side 

92.2 

177.9 

10.3 

Tail 

84.0 

158.6 

13.5 

Starboard 

92.2 

177.9 

10.3 
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T 
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of Boeing 727-200, with no other scatterers present. The received signal is simulated 
by calculating the four beams radiated by the TCAS array. By reciprocity, the 
transmitting and the receiving patterns of the array are identical [22]. The azimuthal 
space is divided into four quadrants, namely 0°-89°, 90°-179°, 180°-269°, and 270 - 
359°. The two beams with the largest amplitude in any given quadrant are known 
from the received signals, and one of these is selected as a reference signal. The other 
signal is subtracted from this reference signal to obtain a value that can then be used 
in the lookup table. Since there are four beam headings, there are four monotomcally 
increasing sections of the lookup table as shown in Figure 2.5. The channel numbers 
of the two largest received signals are also stored. The list of channel numbers and 
the reference signal for all possible <f> is given in Table 2.2. 

Next, the radiation patterns of a TCAS antenna on an actual aircraft are com- 
puted taking into account the effect of various structures on the airframe, like wings, 
tail, engine, etc., which distort the signal transmitted/received by the TCAS array. 
These distortions are in the form of variations in amplitude and phase of the received 
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Table 2.2: Determination of the reference beams for the TCAS under ideal condi- 
tions. 


4 > 

(degrees) 

Highest 

Channel 

Next 

Channel 

Reference 

Beam 

0-44 

0 

1 


45-89 

1 

0 

1 

90-134 

1 

2 


135-179 

2 

1 

2 

180-224 

2 

3 


225-269 

3 

2 

3 

270-314 

3 

0 


315-359 

0 

3 

0 


Table 2.3: Location of intruder from received signals under ideal conditions. 


Reference 
Beam # 

Section of Lookup 
Table searched 

i 

0-89 

2 

90-179 

3 

180-269 

0 

270-359 


signals. Again, the difference signal is computed and stored along with the channel 
numbers in what is defined as the wing table. 

To calculate the bearing error, one starts with a given difference signal and the 
channel numbers from the wing table. From the channel number of the reference 
beam, the quadrant of the intruder can be determined. This is given in Table 2.3. 
That section of the lookup table is then searched for the same value of the difference 
signal as the received signal as depicted in Figure 2.6. The predicted bearing is 
found via linear interpolation when the difference signal lies between two points in 
the lookup table. The bearing error e is defined as 

e = fa - 4>a ( 2 . 6 ) 
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Figure 2.6: Calculation of bearing error. 

where, <f>j and <f> a are the detected (predicted) and actual azimuth angles of the 
intruder, respectively. 

There are two sources of ambiguity/uncertainty that arise in this calculation. 
One type of ambiguity occurs when the difference signal lies outside the range of 
the lookup table. In this case, the two highest signals are a valid combination of 
beams, but the value of the difference signal is not within the boundaries of that 
section of the lookup table. This generates an out of range warning, and the error is 
calculated by setting the difference signal equal to the closest (extreme) value of the 
lookup table. The other source of error arises due to structural scattering, when two 
beams that have opposite bearings have the largest amplitudes in a given quadrant. 
This leads to an illegal combination of channels which leaves the processor with no 
lookup table to search. Hence, one can not estimate the bearing of the intruder. 
The error is then set equal to the previous value of bearing error and a no lookup 
table warning is issued. 
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SPA antenna 


Figure 2.7: Geometry for the spiral phase antenna. Note the antenna alignment. 

2.3 Spiral Phase Antenna System 

As the name implies, the Spiral Phase Antenna (SPA) TCAS operates using the 
phase information contained in the received signal. Note that spiral refers to the 
phase of the received signal, as opposed to the type of antenna. The feasibility of 
using a spiral phase antenna for direction finding applications was first explored 
in [23]. The SPA array consists of four short monopoles arranged on a circle as 
shown in Figure 2.7. Note that this array is offset by 45° in comparison to the CRA 
antenna. As depicted, the positive x-axis points towards the nose of the aircraft and 
the distance d = 2.709" corresponds to A/4 at 1090 MHz. 

The operation of the SPA TCAS is fairly simple. The TCAS equipped aircraft 
transmits an omnidirectional interrogation (query) signal. The transponder on the 
intruder responds by transmitting a reply at 1090 MHz, which is received by each 

element of the TCAS array. 

For modeling purposes on the digital computer, the received signals are simulated 
by transmitting signals at 1090 MHz. In contrast to the CRA antenna, aU monopoles 
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of the SPA array are excited with exactly the same amplitude and phase in the 
airborne antenna code. Sum (E) and difference (A) signals are generated according 


to: 


Ef 

Et 


3 


J2 Ee, and 

i=0 

(E$ 0 — Ee 2 ) + j(Ee 3 


Ee,) 


(2.7a) 

(2.7b) 


where, Ee,,i = 0, 1,2,3 are the complex received signals associated with the four 


elements of the antenna array at some azimuth angle. 

Consider four monopoles of length // 2 on a ground plane, arranged as shown in 
Figure 2.7. The pattern is still described by Equations (2.3) and (2.4). For the SPA 
antenna, the distance between adjacent elements is d, = A/4. Hence the radius a of 
the array is given by 



( 2 . 8 ) 


For the SPA system, I, = A is constant. Hence, from Equations (2.7), one finds 


that 


0 


and 

Et 


j 

A e jkas\r>6cos(4>-4>i) 


i== 0 


4 A cos (— sin # cos j cos (— sin 8 sin <pj 


(2.9) 


2 jA sin ^ sin 0(cos <f> + sin (j>)j — 2 sin sin ^(cos <f> — sin (f>)j (2.10) 


In the lit region of the antenna, due to the size of the antenna array and distance 
to the intruder (who is in the far zone of the array), the amplitudes of all received 
signals are equal; they differ only in phase. The phase angle a is defined by 

a = Mg (|) (2.U) 

In this algorithm, the phase of the sum signal, E% serves as a phase reference. Note 
that a is monotonic in nature and hence can be used to unambiguously determine 
the angle of arrival of the received signal. 
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As is the case with the amplitude system, one must first, generate a lookup 
table. This is done by calculating the radiation patterns of the TCAS antenna on 
the bare fuselage of a Boeing 727-200, located 380" from the nose. The patterns 
are calculated at 10° elevation for a top mounted antenna and —10° for a bottom 
mounted antenna as before. The amplitude and phase of the sum and difference 
beams for a top mounted antenna at the above location are illustrated in Figures 2.8 
and 2.9, respectively. 

It is noted that the pattern origin is at the center of the array and hence the 
phase of the sum beam is constant over all azimuth. The phase of the A beam is 
approximately linear. The monopulse curve (lookup table) is shown in Figure 2.10, 
which also exhibits the monotonic, almost linear variation of phase. This phase curve 
is the calibration curve used in all further calculations involving the SPA TCAS. It 
remains the same even if one is modeling a different airplane. 

Once the lookup tables have been obtained, the radiation patterns are calculated, 
taking into consideration, the scattering and distortion by various structures on the 
airplane. These monopulse curves are distorted due to scattering by the structures 
on the airplane and coupling from other nearby antennas, if present. As in the case 
of the amplitude system, the error curves are calculated for the locations, elevations 
and airframes of interest. 

2.4 Results 

Typical results are provided for the case of a top mounted TCAS on a Boeing 737- 
300 based on the algorithms described above. A scaled line drawing of the principal 
views of the Boeing 737-300 is illustrated in Figure 2.11. As seen, this aircraft has 
two engines mounted below the wings. 

The computer model used for the top mounted antenna is shown in Figure 2.12. 
The engines were not included in the model because of their distance from the 
antenna and the shadowing of the electromagnetic rays by the fuselage. The engines 


18 



f'!"'' 2 8: Ampli ‘ u ^ and phase of SUI " (£) of an SPA antenna. The TCAS 
Elevation™ 10" ^ ° n t<>P ° f * B ° ei " 8 727 ‘ 200 ’ 38 °" from ll >' nose. 
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Figure 2.9: Amplitude and phase 
same location as in Figure 2.8. 


of difference (A) beam of the SPA antenna at the 
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do play an important role in the performance of the bottom mounted antenna. 
The effect the engine has on the TCAS is discussed in Chapter 5. The results 
presented here thus represent the typical errors one might expect on an aircraft 
fuselage when the engines are not visible. It is noted in passing that, normally, 
the contribution from the horizontal stabilizers is often ignored when the radiating 

element is a vertical monopole, as in this study. 

The lookup tables for the CRA and SPA systems were generated on the bare 
fuselage of a Boeing 727-200. The antennas were located 380" from the nose. The 
elevation of the receiver was 10° for the top mounted antennas. Error curves for 
the Boeing 737-300 aircraft are shown in Figures 2.13-2.18 for the TCAS antennas 
located about 289" from the nose (location code 0240). These error curves cover 
elevations —20° to 30° in 10 steps. 

One notices that the error curves for the CRA and SPA systems are surprisingly 
similar, even though the bearing algorithms are completely different. This empha- 
sizes the fact that the bearing error computed here is due to the body of the airplane; 
since it is the same for both systems, the error curves are also similar. The error is 
minimum around 10° elevation. This is because the lookup tables were computed at 
that elevation also. As the elevation goes higher or lower, the error tends to increase. 
This is especially pronounced in the error curves at elevations -10° and -20° where, 
the intruder lies in the shadow region of the TCAS antenna. One may also notice 
that the errors of the SPA system tend to be higher than the CRA system, because 
it is a phase system, and the received phase varies much more rapidly than the am- 
plitude. Finally, due to the way in which the error is calculated, it is seen that the 

error curves of the CRA system are odd-symmetric (e(<£) = -e(-<f>)). 

It is seen from the top of Figure 2.17 and 2.18 that there are discontinuities in the 
error curves of the CRA system at angles where the beams switch, i.e., 0°, ±90 , and 
180°. The reason for this has already been explained briefly. Here, let us attempt 
to elucidate by considering an example. The reader is referred to Figure 2.19, which 
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Figure 2.15: Error curves for CRA (top) a 
TCAS setup same as Figure 2.13. 
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is an enlarged view of the azimuth region 80° to 100°. 
table at 30° are depicted. 


The lookup table and wing 


Case 1: (f> a = 88°. 


Let the actual bearing of the intruder be fa = 88". The channels of the two 
highest received signals are 0 and 1. From Table 2.3, it is determined that the 
intruder lies in the first 0»-90" quadrant. Searching the corresponding section of 

the lookup table, for the value of the received difference signal (9.4 dB), it is found 
<f>d = 82.6°. Therefore, the error is e, = <f, di - ^ o) = -5.4°. 

Case 2: (p a = 89°. 


Let the intruder be at = 89" for the same situation considered above. Now, 
the difference signal is -9.4 dB, and the two highest channels received are 1 and 
2. Then, searching the lookup table to the right side in Figure 2.19, the detected 
position Of the intruder is 97.2" to lead to an error of t2 = ^ = 8 , r . H enc=, 

the magnitude of the discontinuity is: £;-£, = 13.6", and this is easily verified from 
either Figure 2.18 or 2.19. 

Now, it is clear that the jumps in the CRA system occur due to the beam 

switching scheme and the slight change of the beams as a function of elevation. 

By a similar reasoning, it is deduced and verified that there are no discontinuities 

in some error curves, especially close to the elevation at which the lookup table 

was created. Further, the magnitude of the discontinuity increases with increasing 

elevation. This is probably one of the most serious drawbacks of the four element 
CRA system. 

In all cases studied, and those presented here, it is observed that the magnitude 
of the error and the region of its occurrence correlate well with the relative locations 
of the TCAS antenna and the scatterer. This is seen in all figures, where, in the tail 
region, the error is much larger due to scattering by the vertical stabilizer. 

The performance of both systems is similar in the absence of scatterers; however, 
when scatterers are present, the CRA tends to perform better than the SPA TCAS.’ 
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,rom ‘ he 

nose at 30° elevation. 



The implications of this phenomenon will become more apparent in the 
chapters. 


ensuing 
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Chapter 3 

Optimal Location of TCAS 

3.1 Introduction 

One of the important applications of the Bearing Error Transfer Functions (BETF) 
is in the estimation of an optimal location for the top and bottom mounted TCAS 
antennas on aircraft. In this study, three representative commercial airframes are 
chosen, based on the location of engines and the number of planes in service. These 

are: 

Boeing 727-200: Class-3 engine Turbojet, approximately 1300 aircraft in service. 

Boeing 737-300: Class-2 engine Turbojet, about 800 aircraft in service. 

Boeing 747-200: Class-4 engine Turbojet, approximately 180 aircraft in service. 

The criterion used here is based on minimizing the standard deviation of the 
BETF over the angular sectors given below: 

. Top Antenna: Forward sector, -90" to +90" azimuth averaged over -10" to 
-f 20° in elevation, in 5° steps. 

. Bottom Antenna: 360'azimnth, averaged over -20" to 0" elevation, in 5"steps. 

For both top and bottom mounted antennas, the error curves are calculated 
over the whole azimuth with a resolution of 1" at a chosen elevation. The standard 
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deviation /op ( 07 ) and bottom ( 07 .) antennas are defined by 


o-t(6) = 


N ^t(b) + -M|(6) ~ 1 


(3.1) 


where fit(b) is the average value of the BETF averaged over the appropriate range of 
azimuth and elevation angles, that is, 

1 


Pm = 


^h(fc) ^t(b) 

£ £ M 


(3.2) 


N ui') + Mm, ,tr “i 

where e,j is the bearing error and i and j denote the i ,h and j lft elevation and 
azimuth angles, respectively. It is noted that the absolute value of BETF is taken 
in the above equations because based on the definition of the bearing error, BETF 
can be positive or negative. Mt(b) denotes the number of elevation angles for the 
top (bottom) mounted antennas, respectively. Likewise, N t ^ denotes the number of 
azimuth angles for the top (bottom) antennas, respectively. In this study, M t = 7, 
M b = 5, N t = 181, and N b = 361. 


3.2 Boeing 737-300 

The aircraft first considered is the Boeing 737-300. This airplane has two engines 
mounted below the wings as seen from the scaled line drawing of the principal views 
of the Boeing 737-300 in Figure 2 . 11 . The engines do not play a significant role in the 
determination of optimum location for the top mounted antenna; hence, they can 
be ignored in the high frequency model. The best location for the TCAS antenna on 
the bottom, however, is sensitive to the distance from the engine; hence the engines 
must be included in the UTD model. 

The airplane model used for the top antenna has already been introduced in 
Section 2.4. The performance of both the CRA and SPA systems is assessed by 
calculating the standard deviation curves. The standard deviation curve for these 
antennas is shown in Figure 3.1. The detailed statistics of each location may be 
found in Appendix A. The optimal location of the top mounted antenna is at 289" 
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and 229" (location codes 0240 and 0180 respectively) for the CRA and SPA systems, 
respectively. It is seen that the standard deviation curves are fairly constant for both 
systems. This means that one may place the TCAS antenna anywhere from about 
220" to 520" from the nose. The optimal locations are illustrated in Figure 3.2. 

The bottom mounted TCAS antenna is considered next. The computer mod- 
els used for the bottom mounted antennas on the Boeing 737-300 are shown in 
Figures 3.3 and 3.4. It is noted that the aircraft models are different for top and 
bottom mounted antennas and also for the forward and aft region of the airplane. 

Unlike the top mounted antenna case, the engines play an important role in 
the determination of the optimal location because they scatter energy, and hence 
affect the accuracy of the TCAS system in estimating the AOA of the intruder. A 
worst case calculation is made in this chapter with the engines being modeled by 
closed boxes. A flat plate is placed at the front face of the engine to cover the 
inlet. It is important to note that this model of the engine will overestimate the 
specular reflection component of the scattered field. This leads to unnaturally high 
and low errors in the individual error curves. The standard deviation curves given in 
Figure 3.5 confirm this. It is seen that the standard deviation is largest in the region 
from 400" to 600" from the nose, which is also the region occupied by the engines. 
Thus, the optimum location of both systems occurs 95" from the nose, which is 
also the farthest region forward of the engines. This location is shown marked in 
Figure 3.6. It may also be observed that the standard deviation decreases as the 
antenna is moved away from the engines. As in the case of the top mounted antenna, 
the overall performance of the CRA TCAS is better than the SPA. 

3.3 Boeing 747-200 

This aircraft, one of the biggest passenger carriers, poses some unique modeling prob- 
lems. This is due to the presence of a cupola (hump) on top of the forward section 
of the fuselage, which approximately extends up to station 900. (See Figure 3.7). 
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Figure 3.2: Location of optimum position for top mounted CRA and SPA arrays 
Boeing 737-300. 
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Figure 3.4: Computer model for Boeing 737-300 for a bottom mounted antenna 
located aft. Bottom view is seen. 
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Figure 3.5.- Standard deviation for bottom mounted TCAS array: Boeing 737-30< 
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Figure 3.7: Scaled line drawings of the Boeing 747-200 aircraft. 

Thus, a more sophisticated modeling scheme was adopted ineolving theHeUcopte, 
Antenna Radiation Program (HARP) (24). The cross sections of the Boeing - 
fuselage were specified at various locations along the length, as well as the geome . 
of the wings, engines, tail, etc. The code then generates the high frequency compu er 
mode, for the aircraft for an antenna located anywhere on the fuselage. 

For obvious reasons, the mode, for the top of the fuselage for antennas mounte 
on or behind the cupola are different. Two representative models are shown n 
Figures 3.8 and 3.9. The part of the fuselage that is not modeled by the compos, e 

.,Ls » — <•' *• •“ “ 

alL a Z «aa a. a. » — — ‘™“ “ " " 

on the radiation patterns of the TCAS. 


Figure 3.8: Computer model for Boeing 747-200: top mounted antenna (forward). 



Figure 3.9: Computer model for Boeing 747-200: top mounted antenna (aft). 
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The standard deviation curves for the ORA and SPA TOAS systems are shown 
in Figure 3.10. Note that station number 90 is at the nose of the aircraft. The 
corresponding statistical data is given in Tables A. 5 and A. 6. It is observed that the 
standard deviation is very small and similar for both systems. The optimal location 
for both antennas turns out to be at station 0540, where the monopole antennas are 
approximately vertical with respect to the fuselage. The optimum locations for the 
top mounted TCAS are shown in Figure 3.11. Since the standard deviation curve is 
fairly constant, the TCAS antenna may be placed at any convenient location between 
stations 0400 and 0800. This region lies on top of the cupola. If the antenna is moved 
farther back to stations 1380 and 1440, the standard deviation is much larger. This 
is to be expected, and arises due to scattering by the cupola, which now lies directly 
in front of the antenna. 

As in the case for the Boeing 737-300, the engines have to be included when 
simulating the TCAS mounted on the bottom of the fuselage. Two typical computer 
models for the aircraft for bottom mounted antennas are shown in Figures 3.12 
and 3.13. As mentioned earlier, these are the worst case error calculations, in that the 
engines are modeled by closed, perfectly conducting boxes. The standard deviation 
curves as a function of station number are given in Figure 3.14. The optimum 
locations for both systems happen to be at station 0320 and are shown marked on 
a line drawing of the airplane in Figure 3.15. 

The engines do not affect the performance as badly as the Boeing 737-300. The 
reason for this is twofold. First, the engines are farther away from the TCAS antenna, 
and second, the engines are slightly higher than the level of the bottom of fuselage. 
Since elevations in the range of 0° to -20° are considered, the scattering from the 
engines affects the patterns somewhat less. As a result, the standard deviation is 
seen to increase as the antenna approaches the mid section of the fuselage, but the 
rise is not as steep as the Boeing 737-300. Again, the CRA system performs better 
than the SPA system for the locations considered. 
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Figure 3.10: Standard deviation for top mounted antenna: Boeing 747-200. 
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Figure 3.11: Optimum location for CRA and SPA top mounted arrays: Boeing 
747-200. 




Figure 3.12: Computer model for Boeing 747-200: bottom mounted antenna (for 
ward). 
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Figure 3.14: Standard deviation for bottom mounted antenna: Boeing 747-200. 
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Figure 3.15: Optimum location for CRA and SPA bottom mounted TCAS arrays 
on the Boeing 747-200. 
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Figure 3.16: Scaled line drawings of the Boeing 727-200 aircraft. 


3.4 Boeing 727-200 

The Boeing 727-200 airframe was included in this study for its rather unusual place- 
ment of engines two on the fuselage on either side, and a third engine directly 
underneath the vertical stabilizer on top of the fuselage. The scaled line drawings 
for this aircraft are shown in Figure 3.16. The computer model for the top mounted 
antennas is shown in Figure 3.17. The corresponding standard deviation curves are 
given in Figure 3.18. The most striking feature of this error curve is the relative 
flatness of the CRA curve, as opposed to the steep increase in the SPA error. The op- 
timal locations for the CRA and the SPA systems are at 200" and 140", respectively, 
from the nose. These locations are marked on Figure 3.19. 

It is important to note that the space available on the fuselage of a commercial 
airplane to locate a TCAS antenna is very limited due to the presence of many 
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Figure 3.17: Computer model for Boeing 727-200: top mounted antenna. 

other antennas. Thus, if the location of the TCAS antenna were based solely on 
the performance depiction in Figure 3.18, the CRA system has an edge over the 
SPA because its standard deviation is smaller and more constant. Therefore, the 
CRA antenna may be located anywhere between 100" and 580" and still keep the 
standard deviation of bearing error less than about 1.5°. The SPA antenna, in 
contrast, is restricted to be in the far forward region of the aircraft, because its 
standard deviation increases rapidly. This is partly due to the specular reflection in 
the forward region from the plate covering the front face of the center engine. It is 
emphasized that, this is the worst case calculation and a better solution, based on 
the self-consistent multiple scattering matrix formulation is presented in Chapter 5. 

The bottom mounted antenna on the Boeing 727-200 behaves differently. The 
high frequency model used for this computation is shown in Figure 3.20. In contrast 
to the top mounted antenna, the performance of both systems is very similar as 
shown in Figure 3.21. The corresponding optimum locations are shown marked on 



Figure 3.18: Variation of standard deviation with location for top mounted TCAS 
array for the Boeing 727-200. 



Figure 3.19: Location of optimum position for top mounted CRA and SPA antennas 
on the Boeing 727-200. 
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Figure 3.20: Computer model for Boeing 727-200: bottom mounted antenna. 

Figure 3.22. It is seen that the antenna on the bottom may be placed anywhere from 
about 220" to 380" and still have a standard deviation close to the value reached at 
the optimum location. 

3.5 Discussion 

From the above simulations of both systems on various aircraft, some general trends 
are apparent. These are discussed below. 

The error generally increases as one moves the antenna closer to the nose. This 
is very noticeable in Figures 3.1 and 3.10, for both CRA and SPA systems, and for 
the SPA system in Figure 3.18. This is due to the increasing slope of the fuselage 
as the distance between the TCAS antenna and the nose decreases. Recall that the 
calibration table was created for a location of the antenna that is relatively flat; i.e., 
the antenna normal points straight up. When the TCAS array is very close to the 



DISTANCE FROM NOSE ( inch) 

Figure 3.21: Standard deviation curves for bottom mounted antenna on the Boeing 
727-200. 
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Figure 3.22: Location of optimum position for bottom mounted CRA and SP 
antennas on the Boeing 727-200. 





nose the anient,* normal is pointing towards the front of the aircraft. This tdt can 
be as great as 25” as in the case of the Boeing 747-200 for station 0270, as » orpre 

from liable A.5. Thns, there are farger variations in the asimnth pattern than for 

straieht up This leads to a higher error 
locations where the antenna points more s g P 

being competed for locations closer to the nos. than those farther aft. 

U •„ aiso seen that the error is fairly constant for both TCAS systems when^y 

are away from structural scatterers. The effect of the engine is not.cea e 

carves for bottom mounted antennas for the Boeing aircraft, 737-300 an - , 

and for the top mounted antenna on the Boeing 727-200. The role of these engine 

; a m ove the optimum location forward. Note that the effects o, the engine te ds 

. , , ,H h the rising portion of the bottom due 

to decrease as the antenna ,s far forward, in the nsing p 

to the natural blocking offered by the fuselage. As the antenna moves 

! ' se Thus for the bottom antenna, the optimal TCAS location is a 

fuselage increase, lnu , discussed 

“"c«r:i s . v»- - — - - 

— - » 

phase error of about 33.2° /inch of error in distance measuremen 

the CRA system, which operates by comparing the amplitudes o the «ce, e 

* - - — 

much less for errors in distance measuremen , 

M r - . Thus, for a given received signU, the amplitude system is less prone to error 
in calculation as well a. measurement. by sparing the 

Recall that the bearing error for a TCAS location , 

received signals in a lookup table. The lookup table used in t e case ^ 
system has four distinct, monotonic sections, one - „ 

chosen depends on the largest two signals received. The range 
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about 20-24 dB over each section, which is about 9<F wide. Hence, considering the 
worst case, an amplitude error of 1 dB result, i„ an angle of arrival error of 4 5°- 
but, for typical cases, the amplitude errors are much smaller, and hence the error 
tn the AOA. The SPA system also ha, a lookup table which is mouotouic over the 

whole azimuth (slope is unity). Thus, a phase error of 1" !eads to the same amount 
of error in bearing determination. 

This affects the error calculation as follows: Consider the received signal from 
• g-ven scatterer at a certain distance. The procedure, to determine the reflec- 
tion/diffraction points on a plate, geodesic paths on the composite ellipsoid etc, 

are mostly by numerical search. The receiver! vi j 

ine received signal amplitude remains constant, 

whereas, the phase is affected by the distance to tK« a 

uy me distance to the scatterer, as well as the loca- 
tion of the specular point on the scatterer Time l < . . 

ne scatterer. lhus, large phase variations occur due 

to small numerical errors in distance measurement. 

Similarly, when the receiver is shadowed by the fuselage, a numerical search is 
conducted to hnd the location of the effective source on the fuselage, from where 
e ray ,s launched tangentially to reach the receiver. For the same reasons cited 
above, the performance of the SPA system degrades much more rapidly than the 

5y!lem ’ “ lh ' receiw mOT ' S i" ‘bo shadow region, or as the geodesic 

path gets longer. 

Finally, consider an intruder in the forward quadrant of the airplane. Consider 
also, a scatterer, say a vertical stabilizer or a blade antenna to the aft of the TCAS 
For the CRA system, the received signals combine to yield the highest channels as 
and 1. The scattered signal from the blade is very small and it will affect 
the bearmg determination. On the other hand, with the SPA system, there is no 
such beam switching employed, and the transponder signal from the intruder in the 
forward region receives the same weight as a stray scattered signal. Hence, even 
t ough the SPA and CRA systems have comparable performance in the absence of 
scatterer,, the CRA system performs better in the vicinity of scatterer,. 
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Chapter 4 

Effect of Other Antennas 


4.1 Introduction 

In the previous chapters, the operation of the TCAS system in the absence of any 
nearby antennas was considered. The distortion of the radiation patterns by the 
airframe resulted in an incorrect bearing being computed. There, it was assumed 
that the TCAS antenna was far away from (many wavelengths) all other antennas. 
This is not a realistic assumption on a modern commercial aircraft, where one might 
expect to find up to one hundred antennas, all vying for generally the same area on 
the fuselage. When the TCAS is in operation, it induces currents on the antennas 
located in its vicinity. These induced currents then radiate electromagnetic fields 
and distort the radiation pattern of the TCAS array. The magnitude of the induced 
current depends upon many factors, such as distance to the second antenna from the 
TCAS array, size of the antennas, power radiated etc. Moreover, when the frequency 
of operation of the other antenna is close to that of the TCAS, the actual induced 
currents in the TCAS array can be very different from the design. Thus, a smaller 
antenna might also have a significant effect on the operation of TCAS. Hence it is 
necessary to model each antenna in the vicinity of the TCAS and study the changes 
in the behavior of the TCAS as a whole. 

Two blade antennas are considered here because they are found on all com- 
mercial aircraft. Mathematical models are developed to include the effect of the 
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electromagnetic re-radiation from these nearby antennas on the TCAS currents. 
The performance is re-evaluated based on the bearing error curves. It is seen that 
the nearby antennas can (and often do) significantly affect the performance of the 
TCAS. Finally, it is also seen that the effect of nearby antennas on the SPA TCAS 
is much worse than that of the CRA system. 

4.2 Modeling a Nearby Antenna 

Two blade antennas, which are the most common antennas used for VHF and L- 
band communications, are likely to be found in the vicinity of the TCAS antenna 
and are chosen for this study. These are the following: 

1. L-Band blade antenna (Mode S): Dome and Margolin Model AT-741. 

2. VHF communications antenna: Dome and Margolin Model DMC60-1. 

The line drawings of these blade antennas are depicted in Figures 4.1 and 4.2. As 
can be seen, both antennas are thin blades due to aerodynamic considerations. The 
AT-741 is only about A/4 in height; whereas the DMC60-1 is over one wavelength 
tall at the TCAS receiving frequency (1090 MHz). 

In the UTD analysis using the NEWAIR code, one can synthesize a fairly arbi- 
trary variety of radiation patterns by a set of short monopoles and/or slots. However, 
it is necessary to first determine the radiation pattern of the antenna under study. 
Moreover, the above code can not calculate the mutual coupling between various 
antennas mounted on the fuselage; thus, it can not calculate the currents induced on 
the surface of the blade by the TCAS array. Hence, it is necessary to first calculate 
the currents induced on the blades by other means, and next, the fields radiated 
by these currents. Once the fields of the nearby antennas are found, they can be 
superimposed with the radiation pattern of the TCAS antenna to include the effect 
of the blade. 
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The scheme used to accomplish this task is as follows. Since both antennas are 
normally mounted near the centerline of the fuselage, they are within the paraxial 
region. Thus, one can not use the high frequency methods that compute the mutual 
coupling between the two antennas (25], Since the curvature of the fuselage is very 
small along the axis and the antennas are close to each other, one can approximate 
the fuselage by a flat surface. One may then model the arbitrary geometry via 
the moment method (MM). For this case, where the radiation pattern of the blade 
antenna has to be calculated, the Electromagnetic Surface Patch Code (ESP4) is 
used to model the geometry in free space. To this end, the TCAS array is modeled 
by a set of short monopoles. The voltages at each feed point are adjusted to generate 
modal currents of the same magnitude and phase as the currents of the model used 
in the UTD analysis in the NEWAIR3 code. The radiation pattern is then calculated 

via both MM and UTD to verify the accuracy of the model. 

Next, each blade is modeled by a flat plate. Based on measurements conducted 
at MIT Lincoln Laboratory, it was determined that the resistive termination (load) 
placed at the terminals of the blade had very little effect on the radiation patterns 
of the TCAS antenna (26]. From this, it may be concluded that the nearby blade 
antenna behaves much more like a scatterer than an antenna. Hence, the blade 
could be modeled in each case as a scatterer. Once the currents induced on the 
blade surface are calculated for a particular separation and elevation, they are stored 
for later use. Since the steps used in the analysis of L-band and VHF blades differ 
widely, they are discussed separately in the following sections. 

4.3 L-band AT-741 Blade 

Keeping in mind that only short monopoles, normal to the fuselage at the point of 
attachment can be used in the NEWAIR3 code, one first replaces the blade by a set 
of four wires strategically placed as shown in Figure 4.3. This is possible because 
the dominant induced currents on the blade are predominantly vertical, such that 


56 


AT741 wire model 
Dimensions in Inches 


[-* o. 9*4—0. 9*4— o. 9^ 


Figure 4.3: Wire model of L-band blade used in the moment method analysis. 

the horizontal component can be neglected. This is demonstrated in Figure 4.4, 
where the patterns obtained via the moment method for the blade model and the 
four wire model are shown for a blade located 4 feet (48") aft of the CRA antenna 
array. The blade/wire model was illuminated by the beam pointing towards the 
blade (beam 2), and the radiation pattern was computed at 10° elevation. One 
should observe that the two radiation patterns are in good agreement. Since the 
monopoles used in the UTD analysis have to be under A/4 long, the two wires in the 
middle, which were slightly longer than A/4 were trimmed to A/4 and the induced 
currents and radiation pattern were computed again to verify that they did not 
change appreciably. The final wire model for the AT-741 blade used m the UTD 
analysis is shown in Figure 4.5. 

The final step in this synthesis process is the inclusion of the radiation pattern 
of the blade in radiation patterns for the TCAS array. This can be accomplished via 
two methods. One can either recompute each element pattern in the presence of the 
blade antenna and process the received signals, or, compute the blade scattered fields 
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Top 


view 




Front view 


Wire model for AT741 

Figure 4.5: Wire model used in NEWAIR3 code for the AT741 antenna. 

by illuminating it by the four CRA beams or by the sum and difference beams of the 
SPA TCAS in turn. These can be then directly added to the corresponding received 
signals. (Of course, one must be careful about the phase centers of the TCAS and 
blade antennas). The latter approach is adopted here, due to its flexibility and the 
fewer computations required. 

In case of the CRA array, there are four beams each pointing in a different 
direction. Hence the illumination of the blade is different for each beam, and so are 
the induced currents. Thus, the CRA beams are created via the moment method by 
weighting the monopoles correctly, and the induced currents due to each beam on 
the wires are calculated and stored for use with the NEWAIR3 code. Similarly, for 
the SPA array, the sum and difference beams are created separately and the induced 
currents on the wires due to this type of illumination are also calculated and stored. 

The curvature of the fuselage has only a small effect on the patterns of the 
blade antenna, because both antennas are along the axis of the fuselage, where the 
curvature close to zero. Hence, the approximation of the fuselage by a ground plane 
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Figure 4.6: Radiation pattern of wire model for the L-band blade mounted on the 
fuselage of Boeing 727-200, calculated via UTD. 

is valid. To demonstrate this, the wire model of Figure 4.5 was mounted on top of a 
Boeing 727-200 fuselage, 380" from the nose. The scattered fields of the blade model 
for the CRA beam pointing towards the blade (beam 2) is shown in Figure 4.6 at 
an elevation of 10°. Note that, even though the levels of the two field patterns differ 
by about 34 dB in the moment method computation and those obtained via UTD, 
the relative levels of radiation from the blade antenna with respect to the peak level 
of the TCAS beam is the same (about 17 dB). 

4.4 VHF DMC60-1 Blade 

The VHF communication blade is one of the most common antennas found on com- 
mercial aircraft. It is even included in the line drawings of the airplanes given in 
Figures 2.11, 3.7, and 3.16. It is a general purpose VHF communication antenna op- 
erating in the 118-136 MHz range. Modeling this antenna for use in the NEWAIR3 
code presents some unique problems, most of which are due to the size and shape of 
the antenna. As in the case of the L-band blade, the plate model obtained by the 
moment method is replaced by a set of inclined wires shown in Figure 4.7. These 
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DMC 60 Blade 
Wire Model 



Figure 4.7: Model for the VHF blade using inclined wires. 

wires are inclined at about 70° with respect to the ground plane. The radiation 
pattern for the flat plate model and the inclined wire models agree well as shown 
in Figure 4.8. The radiation patterns do not agree well when the wires are vertical. 
This indicates that the induced currents on the blade have significant horizontal as 
well as vertical components. The inclined wire model can not be directly used in 
the UTD procedure. The monopole antennas used in the NEWAIR3 software must 
meet the following criteria: 

a. the monopoles must be perpendicular to the fuselage, and 

b. the length of the monopoles must be less than A/4. 

The inclined wires are about one wavelength long at 1090 MHz, and hence, simply 
replacing the wires by short A/4 monopoles will not give acceptable results. Hence 
other methods of synthesis must be explored, which meet the above constraints. 

The objective here is to synthesize an arbitrary radiation pattern denoted by E™ 
using an arbitrary planar array of N vertical monopoles on a ground plane. The 
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Figure 4.9: An arbitrary planar array of N monopoles. 

field E e is the scattered field by the blade, computed in this case via the moment 
method. 

The geometry of the planar array is shown in Figure 4.9. Let the length of the 

i th monopole be defined by l, and its position by and <f>,. The excitation of this 
monopole is given by 


/, = A,e }ai 

where A, is the amplitude, and a, is the phase of the excitation. 

The radiated field of this monopole array in the far-field is given by 


(4.1) 


£7(M) = i; 

i=1 


jvh 

cos 

^ cos ffj 

-cosf 

2tt 


sin 6 



e jkpi sin 6 cos (<£-«£, ) 


(4.2a) 


and 


®*( s .^) = 0 (4.2b) 

at any receiver location (8,*). The snperscript denotes that the field is due to 
the monopole (rare) elements. Now, let E? he given at M receiver positions, 
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N 


where 


£r(^) = £ 7 * z ' j ’ 

»=i 

COS cos 6 j ^ — cos 


7 _ JV 

z ' j - 2tt 


j = 1,2, ...,M 


jfcpi sin 0j c:os( <£j-(£i) ^ 


sin 6j 


(4.3) 


(4.4) 


For M > N, as is normally the case, this denotes an over 
equations. In matrix notation, one finds that 

[ejh = iz][i] 


-determined system of 


(4.5) 


where 


[E?] 
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Zn Z12 
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Z t /v 
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[ Zm\ Z]\I2 9 * * 


n 


(4.7) 


and 


[i] = 


1 1 

h 

[ hi 


(4.8) 


I„ the above equations, [I] is the unknown excitation vector to be solved for. 
This is derived using techniques from linear algebra such that 

(4.9) 


[I] = [([Z] T [Z])] _I [Zf[E?]- 


The solution for [I) thus obtained via (4.9) is the least squares solution 127], Tins 
minimizes the squared difference between the fields due to the blade E? and due to 
the set of wire. Bf. It is noted that the above analysis is fairly general, and can 
be easily modified to account for radiation from an, arbitrary set of antennas, he 
radial and axial slots, and/or a combination of wires and slots. 
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Top view 
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Side view Front view 

Wire model for DMC60 

Figure 4.10: Wire model of VHF blade for NEWAIR code. 

A model that was both accurate and easy to use in the NEWAIR3 code was 
then derived by varying the parameters of the monopole array such as length and 
spacing. It was observed that much better results (smaller number of monopole 
elements, lesser error) were obtained when the above procedure was carried out 
around the phase center of the blade antenna. This is because, when the scattering 
patterns of an antenna are referred to its phase center, variations in the phase are 
minimized. This becomes especially important when the separation between the two 
antennas is large, and the phase varies rapidly. The final model for the DMC60-1 
consisted of six A/4 monopoles arranged as shown in Figure 4.10. 

The model was verified using the NEWAIR3 code by computing the radiation 
pattern of the monopoles located 48" aft of the CRA TCAS antenna, on the fuselage 
of a Boeing 727*200, at 10° elevation. The radiation pattern of this array when 
illuminated by the beam pointing toward the tail (beam 2) is shown in Figure 4.11. 
The agreement is very good. Again, the difference in levels is relative, due to different 
normalizations in the two methods; i.e., the moment method and UTD. 
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Figure 4.11: Radiation patterns of moment method and UTD models of VHF blade 
mounted on fuselage of Boeing 727-200. 

It was found that these two models agree well only for a small range of elevation 
angles around which the synthesis was carried out. This is because, the flat plate 
model has horizontal components of induced currents, which are not accounted for 
by the monopole model. Hence, the synthesized monopole model can not be used 
for more than about 10° around the elevation at which it was created. Therefore, 
unlike the model for the L-band blade, it is necessary to calculate the induced wire 
currents at several elevation angles. As in the case of the L-band (AT741) antenna, 
the induced currents were calculated for each beam of the CRA array as well as the 
sum and difference beams of the SPA antenna. 

4.5 Simulation and Results 

Once the induced currents have been obtained, the monopole models for the L- 
band and VHF blades are placed on the fuselage of the aircraft under consideration. 
Since the total received fields are found using superposition, the TCAS array pattern 
is included as a separate term. The radiation patterns of these simulated blade 
antennas were then calculated in the presence of all scatterers on the airframe and 
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for the various illuminations at different distances from the TCAS array. Finally, 
the corresponding radiation patterns were superimposed, to obtain the patterns for 
the TCAS antenna when it is operating in the presence of the blade. 

Error curves corresponding to a location 289" from the nose on a Boeing 737- 
300, for elevations -10° to 20° were already introduced in Figures 2.13-2.16. In 
these curves, the effect of the vertical stabilizer is apparent in the rear quadrant, but 
overall, the errors are relatively small. Figures 4.12-4.15 depict the same antenna 
location, except that the TCAS is now operating in the presence of the DMC60 
blade antenna. The blade antenna was 48" aft of the TCAS antenna. Similar error 
curves for the TCAS antenna operating in the presence of the AT741 blade are given 
in Figures 4.16-4.19. 

It is seen that the blade antenna has a significant influence on the performance 
of the two TCAS systems. The most obvious effect is the oscillatory behavior of 
the SPA error curve over [ — 180:180] range, as opposed to oscillations only over the 
[-180: -90, 90 : 180] range for the CRA TCAS. This difference is due the way in 
which the received signal from each element is processed by the CRA system. Recall 
that the receiving antennas in this simulation are modeled as radiators. Thus, to 
estimate the bearing of the intruder, the CRA TCAS transmits four beams in the 
assigned directions. The blade, which is between the TCAS and tail, will not affect 
beam 0 except in the backlobe region. Beams 1 and 3, pointing to the sides, are also 
only slightly affected. Only beam 2, which illuminates the blade directly, is severely 
distorted. This is shown in Figure 4.20, where the distortion of the TCAS pattern 
by the VHF blade is apparent. 

Although the blade re-radiates whenever the TCAS is in operation, the radiation 
is strongest only when illuminated by beam 2. Hence, the two highest received signals 
are still via beams 0 and 1 or 0 and 3 in the forward half. The radiation by the 
blade in the -90° to 90° azimuth sector is not strong enough to change the two 
highest signals. Therefore, the distortion of the error curve by the blade is minimal 
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Figure 4.12: Error curves for CRA (top) and SPA (bottom) TCAS in the presence 
of DMC60 blade antenna 48” aft. The TCAS antenna was mounted 289” from the 
nose on Boeing 737-300. Elevation = —10°. 
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Figure 4.14: Error c 
of DMC60 blade, at 
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Figure 4.16: Error curves [or CRA (top) and SPA (bottom) TCAS in the presence 
of AT741 blade antenna 48" aft. The TCAS antenna was mounted 289 from the 
nose on Boeing 737-300. Elevation 10 . 
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Figure 4.17: Error curves for CRA (top) and SPA (bottom) TCAS in the presence 
of AT741 blade, at 0° elevation. Aircraft setup same as Figure 4.16. 
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Figure 4.19: Error curves for CRA (top) and SPA (bottom) TCAS in the presence 
of AT741 blade, at 20° elevation. Aircraft setup same as Figure 4.16. 
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Figure 4.20: Distortion of CRA beam 2 due to re-radiation by the DMC60-1 blade 
antenna. 

in the forward half is small. Of course, the error introduced by the blade in the rear 
quadrant is larger, because the blade is now illuminated by beam 2, which is the 
strongest beam in the rear quadrant. On the other hand, the SPA antenna always 
illuminates the blade via sum and difference beams, which are both omnidirectional, 
and differ in amplitude by only 3 dB. Thus, the effect of the blade is visible over the 
whole azimuth region as shown. 

The behavior of the TCAS in elevation is quite similar for both systems. The 
radiation from the blade drops off approximately as sin# and hence, their error 
magnitude decreases slightly. It is interesting to note that the error magnitude is 
comparable for both the AT741 and DMC60-1 blade antennas. 

Next, the behavior of the two TCAS systems was simulated at various separations 
from the blade antennas. For this study, the TCAS antenna was held fixed at 289 
from the nose on the Boeing 737-300. The blade antenna was placed at increasing 
distances aft of the TCAS antenna, starting at 24" (2.2A) from the TCAS up to 168" 
(15.5A). Error curves were computed via UTD for elevations ranging from 0° to 20 
in 5° steps and averaged to find the standard deviation of the bearing error. The 
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Boeing 737-300. 

results „ { this high frequency computation .re given for the AT24, in Figure 4,, 
Rnd the DMC60-1 in Figure 4.22. 

To verify these high frequency results, error curves were calculate via an i 
pendent moment method computation. For this case, theh.ade antenna was modeled 
Talat piate and the antennas hy thin wires. The curvature of the fusei.ge ^ 
ignored and the fuselage itself was then replaced by an mSmte gr P 

using image theory, the ground plane was removed and the w.re monop ks an 
b ,ad antennas were augmented with their images. The geometry o, the mom 

method mode, for the OT 41 — 

in Figure 4.23. The model for the DMC60-1 is similar and » shown g ^ 
Note that the latter Bgnre depicts CRA array. Radiation patterns were then mm 
ted for each TCAS antenna in the presence of each blade. By processmg the, 
patterns as already described, the standard deviation of the hearing error was 
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Figure 4.23: Moment method model for the AT741 antenna 48" aft of the SPA TCAS 
antenna. 

puted. It is noted here that, due to the modeling procedure, only positive elevations 
are valid. Hence the standard deviation was computed only for elevations in 0° to 
20° range. The process was repeated at various elevations and separations to obtain 
the standard deviation of the bearing error for each location of the blade antenna. 
The standard deviation curves for the AT741 antenna are given in Figure 4.25 and 
the DMC60-1 blade in Figure 4.26. 

The results are very interesting. Firstly, it is seen that the UTD and MM models 
agree. The minor differences in the curves for the high and low frequency models 
can be attributed to the differences in the models used — like flat plates versus 
monopoles for the antennas, flat versus the curved ground plane etc. Secondly, 
the SPA error monotonically decreases as the distance is increased from the TCAS 
antenna. In case of the CRA system, it is small to start with, increases in a region 
about 4 to 6 wavelengths behind the TCAS, and then decreases again. Finally, it is 
seen that, in the case of the UTD models on the Boeing 737-300, given in Figures 4.21 
and 4.22, that the error approaches that of the TCAS antenna alone (dashed line) 
as 6 increases. It is noted that the bearing error due to an isolated TCAS is close 
to zero. The error is identically zero for the intruder at the same elevation as the 
lookup table, and small for other elevations. 
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Figure 4.24: Moment method model for the DMC60-1 antenna 48 aft of the CRA 
TCAS antenna. 



Figure 4.25: Moment method computation to verify the results of Figure 4.21. Fuse 
lage was replaced by a flat ground plane and blade modeled as a flat plate. 



• t verify the results of Figure 4.22. Fuse- 
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TCAS beam is not affected as much “ ^ P “ ^ ^ tbe DMC 60-1 is somewhat 
It is also seen that the error of the llie AT741 is about 3“ 

smaller than the ATT41. This is ^ due l0 the magnitude of the 

ta„; where., the DMC60-! is abou 2 ' ^ ^ ele ctrtcnl sise of 

lanced currents which depend upon the sep 


the blade antenna. 
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Chapter 5 

Effect of Engines 


5.1 Introduction 

Engines play an important role in the performance of the TCAS. The effect of the 
engines is especially noticeable when an engine inlet is directly illuminated by the 
TCAS antenna. To better understand the true scattering of the engine cavity, the in- 
let is modeled by an open-ended circular waveguide cavity in free space. It is assumed 
to be excited by an electromagnetic plane wave. The waveguide is semi-infinite and 
has an interior termination which may be modeled as a flat plate (short circuit), a 
dielectric material, or a flat blade disk optionally with a conical or hemispherical 
hub. Some of these geometries are depicted in Figure 5.1. 

In the analysis presented below, primary emphasis is placed on the circular waveg- 
uide, because it approximates the structure of the engine most closely. Note that 
the terms, engine and inlet, are used interchangeably. In all cases studied here, 
the termination of the circular waveguide is assumed to be a PEC flat plate. The 
aircraft model used in this study is the Boeing 727-200, and the engine on top of 
the fuselage is of interest. Results are presented for both CRA and SPA systems for 
selected elevations in the lit region in the forward half of the airplane. 

It is shown that the model of the engine plays a more prominent role on the 
performance of the SPA system than the CRA. Standard deviation curves to assess 
the systems are also given. These curves clearly demonstrate that the spiral phase 
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Figure 5.1: Some possible waveguide and termination geometries for modeling the 
aircraft engines. 
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TCAS is more susceptible to adverse engine effects than the amplitude system. This 
is especially true when the engine inlet is in the lit region of the TCAS antenna. 
It is seen that the overall errors are slightly smaller, and more uniform. Finally, 
it is noted that the performance in the forward half may be slightly improved if 
the TCAS antenna is closer to the nose, thus utilizing the natural curvature of the 
fuselage to cause some blockage of the engine inlet illumination. 

5.2 Scattering by a Waveguide 

The inlet analysis employs a combination of high frequency techniques such as the 
uniform geometrical theory of diffraction and its modifications based on the equiv- 
alent current method (ECM) as well as the physical theory of diffraction (PTD). 
These are used in conjunction with the self-consistent multiple scattering method 
(MSM) [18]. This method is presented in more detail in Appendix C. The analysis 
is fairly general and has been successfully applied to circular as well as rectangular 
waveguides. 

In the hybrid asymptotic modal method, the self-consistent scattering matrix 
formulation is used to describe the fields inside the waveguide cavity. This includes 
the multiple interactions between the open end and termination. High frequency 
techniques are used to obtain the modal scattering matrices that describe the cou- 
pling between the cavity modes and the external region via the open end, and the 
interior modal reflection from the open end. The modal scattering matrix describing 
the reflection from the waveguide cavity termination may be found in closed form 
for the flat plate or dielectric termination. Physical optics is used to approximately 
compute the reflection matrix from disk/blade/hub terminations. The external scat- 
tering by the rim at the open end is found by replacing the rim with equivalent edge 
currents. These edge currents are found approximately by incorporating the UTD 
diffraction coefficients. First order as well as second order diffraction effects can be 
included in the computation. 
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5.3 Engine Model 


The engine is modeled as a terminated circular waveguide cavity. The scattered 
fields of the terminated circular waveguide are calculated using the MSM formu- 
lation, introduced above. The computation is facilitated by a modified version of 
software titled NEAR_CIRC, which has been developed at the Ohio State University 
ElectroScience Laboratory [20, 19]. This program computes the bistatic scattered 
fields from a rectangular or circular semi-infinite waveguide terminated inside. The 
scattered fields are output as: Em, Eo<f>) 2£<^, andi?^, where the first subscript refers 
to the polarization of the incident plane wave and the second subscript to that of 
the scattered field. It is assumed that the amplitude of the incident plane wave is 
unity. 

Figure 5.2 depicts the geometry of the engine mounted on an airplane. The 
aircraft coordinate system, the TCAS coordinate system as well as the engine’s 
local coordinate system are shown. The coordinate system of the TCAS is the same 
as the pattern coordinate system for azimuth cuts. 

The scattered fields of the engine (inlet) are incorporated into the patterns of 
the TCAS antenna via superposition. This must be done carefully because of the 
various coordinate systems involved. The scattered fields due to the engine are 


Eh 

K,.\ 


iQtwm 


E, 

E: 


Oe 


4>e 


e ~jkff ac 


(5.1) 


Here, Eg t , E^ t are the scattered fields of the engine in the TCAS coordinate system, 
Ek yEfa are the scattered fields in the far zone of the inlet, with the phase referred 
to the center of the engine. [Qt] denotes the transformation matrix from rectangular 
to spherical coordinates with reference to the TCAS coordinate system, such that 


m = 


cos 9 cos <j> cos 8 sin <f> — sin 9 
— sin <f> cos <f> 0 


(5.2) 
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Figure 5.2: Geometry for modeling the engines. 

The spherical to rectangular transformation with respect to the engine center is 
given by 

cos 9' cos (j>' — sin <f>' 

= cos 9' sin <f>' cos <f>' . (5-3) 

— sin 9' 0 

Finally, is a rotation matrix, which for the top mounted TCAS is simply 

given by 

'001 

1 0 0 | . (5.4) 

0 1 0 


[Rte] = 

From Figure 5.1, one finds that 


T : ae — ^ae^ae ® 

e = xx e + yy e + zz e , 
d = xx a T yy a + zz a 


and 


(5.5) 

(5.6) 

(5.7) 


referring to the coordinates of the center of the inlet and the TCAS antenna, re- 
spectively. The incident direction of the plane wave at the inlet is given 
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Finally, note that, for a given pattern direction (0,<£), the corresponding direction 
of the scattered field in the inlet coordinate system is given by </>„), where 


0 S = cos 1 (sin 6 cos <£), and 



(5.10) 

(5.11) 


Eq v and E^ c are the components of the scattered field of the engine (referred to 
engine center), suitably modified to account for the excitation by the TCAS antenna. 
This may be written as 


El 


E. 


4>e 


Eoo Eo<t> 
E<t>e E<t>4> 


Ah 

A 1 


(5.12) 


Here, A l 0e and are the scaled incident amplitudes of the plane wave from 
the TCAS antenna. Eoo, E^ } and E^o are the scattered fields from the inlet, 
with the first subscript indicating the polarization of the incident plane wave and 
the second subscript that of the scattered field. 

Recall that the excitation of the engine is via sum and difference beams for 
the SPA system and the four beams for the CRA TCAS which were derived in 
Sections 2.2 and 2.3. To get the correct amplitude and phase of the incident plane 
wave, one first computes the fields due to the TCAS antenna, and converts them to 
the coordinate system of the inlet. Hence, one finds that 


A 1 

^6c 

A 9 

™4>c 


M \-1 


= ([Qmc)m) 


El 

El 


(5.13) 


where, El,El are the components of the radiated field from the TCAS at the 
engine location. 
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5.4 TCAS Performance 


The performance of the TCAS is re-evaluated in the presence of the engine using 
the above models. This may be accomplished via two equivalent methods. One 
can recompute each element pattern in the presence of the inlet and process the 
received signals as before; alternately, the inlet scattered fields can be computed by 
illuminating it by the four CRA beams or the sum and difference beams of the SPA 
system in turn. These can then be directly superimposed on the receiver beams 
computed for the same geometry without the engine inlet. The latter approach is 

adopted here for its flexibility and lesser computations. 

In case of CRA array, there are four beams, each pointing in a different direction. 
Hence, the illumination of the inlet is different for each beam, and so are the scattered 
fields. Thus, the CRA or SPA beams are created by weighting the monopoles in turn 
to compute the field incident on the inlet. The scattered fields are obtained via the 
MSM method and stored. These fields are then superimposed to obtain the pattern 
of the TCAS antenna in the presence of the engine. Note that, the received signal at 
the TCAS antenna is simulated by computing the transmitted field, which because 
of reciprocity are identical. 

The engine on top of the Boeing 727-200 was modeled via this procedure to check 
the validity of the results. The model of the engine consisted of a circular waveguide, 
terminated by a PEC plate. The radius of the waveguide is 2.89A. The depth was 
2.77A, which is approximately the depth of the blades from the inlet. The TCAS 
antenna was placed successively at 60" intervals at distances 100" up to 400" from 
the nose on top of the fuselage. 

The scattering from the center engine of the Boeing 727-200 is shown in Fignre 5.3 
for 0” and 10° elevations. The incident field at the engine was from beam 2 of a CRA 
TCAS, located 340" from the nose. It is noticed that the scattering from the engine 
is more nniform, due to the large nnmber of modes inside the waveguide structure. 
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Figure 5.3: Scattered fields from a terminated circular inlet at 0° (top) and 10° 
(bottom). The inlet was excited by beam 2 of the CRA TCAS placed 340" from 
nose on a Boeing 727-200. 
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Error curves presented in Figures 5.4 and 5.5 depict, the effect of the engine 
model on the performance of the two TCAS systems at 0° for the engine modeled 
by the circular waveguide and closed rectangular PEC box, respectively. The TCAS 
antenna was placed at location 0340 (same distance in inches from the nose) on top 
of the Boeing 727-200. Similarly, Figures 5.6 and 5.7 and Figures 5.8 and 5.9 depict 
the error curves for the same location at 10° and 20° elevations, respectively. 

It is seen that from these results that the error curves are comparable, especially 
at lower elevations. This is due to the large backscatter from the waveguide ter- 
mination. It is also noticed that the bearing error due to the waveguide model is 
more uniform; i.e., one does not see the specular, localized regions of higher errors 
as in the flat plate model. In the case of the plate model, the highest error drops 
quite rapidly. This is not the case with the waveguide model, whose error tends to 
be more uniform. This is expected because, the waveguide has many propagating 
higher order modes, which give rise to higher side lobe levels. Hence, the distortion 
of the beams is more uniform, and the error tends to be smaller. 

To study if the performance of the SPA system could be enhanced by moving it 
forward, standard deviation curves were generated for each TCAS location spanning 
100" to 400" at 60" intervals along the top of the Boeing 727-200 fuselage. The results 
are depicted in Figure 5.10. It must be noted that, due to modeling constraints, 
only positive elevation are valid. Thus, the standard deviation was calculated in the 
forward sector, -90° to 90° for elevations 0° to 20° only. The figure also shows the 
tilt angle of the antenna normal towards the nose of the airplane. It is seen that the 
minimum error occurs around station 160, where the tilt is about 7.5°. This is close 
to the optimal location at 140, obtained with the plate model for the top engine. 
Comparing Figures 5.10 and 3.18, it is noticed that the overall errors are less for the 
TCAS systems when the top engine is modeled by circular waveguide. The general 
trends are still the same the standard deviation of the error decreases and then 
increases as one increases the distance from the nose. The standard deviation of 


91 
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Figure 5.6: Error curves for CRA (top) and SPA (bottom) TCAS at 10° elevation. 
Top engine of Boeing 727-200 modeled by terminated circular waveguide. 

















bearing error 



Figure 5.9: 
Top engine 


. curves for CRA (top) and SPA (button,) TCAS at 20” elevation, 
oeing 727-200 modeled by closed rectangular box. 


97 



Figure 5.10: Standard deviation curves for Boeing 727-200, with 
by a terminated circular waveguide. 


top engine modeled 


the error i, practically the same tor both cases of engine models for location 0100. 
This was found to be due to the placement of the TCAS. As one moves the TCAS 
forward, the top engine illumination gets blocked by the fuselage. Thus, the incident 
field at the inlet is smaller, because the creeping rays get rapidly attenuated by the 
fuselage. At location 0100, however, the attenuation is so great that the scattering 
from the inlet or the plate modeling the top engine is insignificant. Hence, it is 
concluded that the effect of the engine on the SPA TCAS can be overcome to a 
certain extent by moving the antenna forward. 


98 




Chapter 6 
Conclusion 


6.1 Introduction 

The operation of CRA and SPA TCAS systems has been studied in terms of system 
performance based on electromagnetic effects. Each system has its own unique 
features, and shortcomings. In this chapter, these results are summarized. Possible 
methods of improving the performance are suggested to overcome some of these 
deficiencies. The error increases very close to the nose, due to the forward tilting 
of the TCAS antenna. The error also increases as one moves the antenna aft, i.e. 
towards the stabilizers. 

It is interesting to find that the performance of both CRA and SPA systems 
is similar in the absence of nearby scatterers. The CRA however, performs much 
better in the vicinity of scatterers. Even in the absence of nearby scatterers, the 
standard deviation of the bearing error is consistently smaller for the CRA TCAS 
than the SPA system. A technique that could improve the performance of the CRA 
TCAS is suggested. 

Lastly, the results obtained so far are extrapolated to a couple of other popular 
commercial airplanes, namely the MD-87 and the Boeing 757. It is expected that 
the TCAS would perform well if installed in the forward region of the fuselage based 
on the results generated in this report. 
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6.2 Elevation Effects 


It is now well known that the performance of the TCAS system depends upon the 
elevation of the intruder. In all cases studied, the errors are smallest at 10°, which 
is the elevation at which the lookup table is created. The errors increase as the 
elevation changes. It was also seen in Chapter 2 that the CRA system exhibits 
discontinuities in the bearing error curves at elevations other than the one at which 
the lookup table was created. When the TCAS is in operation, the bearing of 
an intruder is continuously monitored and updated. Some of the signal processing 
involves the rate of change of the intruder’s bearing. When the intruder passes from 
one quadrant to another, the jump in the detected bearing leads to discontinuities 
in the bearing rate signal. Unless this is detected and corrected immediately, the 
TCAS pilot would be given wrong information about the intruder. The spiral phase 
TCAS does not have these discontinuities in the bearing rate curve, because there is 
only one monotonic section of the lookup table, and the phase of the received signal 
always lies in the [—180° : 180°] range. 

The CRA system performance can be easily improved by using elevation depen- 
dent lookup tables. This presumes knowledge of the elevation of the intruder, which 
is available in the Mode-S transponder return. On the other hand, this entails more 
involved signal processing and storage requirements for (possibly) a large number 
of calibration curves, which should not be that big a problem for modern digital 
systems. 

Lastly, it is noticed in general that the performance of the TCAS antenna gets 
worse as the intruder gets deeper into the shadow region of the TCAS (below the 
horizon, for a top mounted antenna) due to the blockage by the fuselage. The 
amount of error induced by the fuselage is fairly constant for the antenna in the 
mid-section of the fuselage, and gets smaller as the antenna is moved forward. 
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6.3 Nearby Antenna Effects 


It was seen that the re-radiation from a nearby scatterer, such as another blade 
antenna, could affect the performance. Two typical blade antennas were considered, 
the AT741 and DMC60-1. Models based on the moment method were derived for 
these blade antennas. These models were then used along with the models for the 
TCAS antenna to study the performance of the TCAS system in the presence of these 
blades. It has also been demonstrated via moment method and UTD simulations 
that the blade antenna patterns are affected only slightly by the curvature of the 
fuselage. Based on these studies, it was found that the SPA system is severely 
affected by the blades. The CRA TCAS does not suffer from this problem to the 
same degree, due to its lobe switching technique. Consequently, it was established 
that the nearby antennas could be much closer to the CRA antenna than the SPA. 
It was also discovered that the nearby antenna need not be physically large to affect 
the bearing accuracy of the TCAS. The AT741 L-band blade, which is physically 
smaller than the DMC60-1 VHF antenna, illustrated this point. 

Thus, for acceptable performance from the SPA antenna, one solution is to move 
it many wavelengths away from the scatterers. This would help two ways; firstly, 
due the increased separation, the re-radiation will be weaker from the other antenna. 
Secondly, the natural curvature of the fuselage would tend to block or attenuate the 
fields reaching the blade antenna. 

The overall performance of both TCAS systems will be severely compromised 
by placing the TCAS aft of the blade antenna. The only solution in this case is to 
either switch the antenna positions (not always possible) or to increase the separation 
between the two antennas. 
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6.4 Engine Effects 

In some of the cases studied, notably the bottom mounted antennas on the Boeing 
737-300, and Boeing 747-200, it was found that for good performance of the TCAS 
antenna, it must be as far forward as practically possible. This is especially true for 
the SPA system, which has been proven to be more susceptible to scatterers around 
it. Since the SPA system does not discriminate as well as the CRA, it was found 
that it is much more difficult to find an acceptable performance location for the SPA 

than the CRA. 

The engine plays a vital role in the performance of top mounted TCAS antenna in 
the case of the Boeing 727-200. It is once again seen that the SPA system is affected 
worse than the CRA TCAS. It was found that the adverse effect of the engine on the 
SPA system could be overcome to a certain extent by moving the antenna forward 
so that the fuselage curvature tends to block the engine illumination. It was also 
found that the CRA system performs much better in the presence of engines and 

other scatterers. 

6.5 Other Factors 

The treatment of the TCAS system has been solely from the viewpoint of the antenna 
designer. It is noted however, that this is last stage of the entire TCAS and and a 
fair amount of signal processing is necessary to acquire and track the other aircraft 
in the vicinity. 

A factor that causes further deterioration of the performance occurs due to 
changes in the characteristics of the antenna, cables, etc., due to temperature and 
other variations. These changes affect both systems; however, it is suspected that 
the SPA system will be more sensitive to phase shifts and hence probably more 
difficult to calibrate and maintain. 
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The accuracy and stability of the phase comparison system may be enhanced by 
utilizing the amplitude information available at the antenna terminals. One could 
combine the antenna outputs with passive circuitry to yield sum and difference 
signals, and process these signals as a conventional amplitude system as suggested 
in [4]. 


6.6 Concluding Remarks 

It has been seen that all these studies generally point to locating the TCAS antenna 
towards the forward region of the fuselage. For the top mounted antenna, it is best 
to locate the antenna where the fuselage is curving down so that good coverage 
can be obtained for lower elevation angles and lessen the effects of the top-mounted 
engines, blade antennas and/or vertical stabilizer. For the bottom mounted antenna 
however, it is best to locate it on a flat section, because full azimuthal coverage is 
desirable. Care must be exercised to mount the antenna away from engines and 
other scatterers. 

Finally, based on these findings, some remarks can be made about other aircraft. 
These are listed below: 

6.6.1 McDonnell Douglas MD-87 

The McDonnell Douglas MD-80 series airplanes, developed from the DC-9 are pop- 
ular as a short /medium range commercial transport. Other aircraft in the series are 
similar, except for the fuselage length. The most recent addition to this family is 
the MD-90 airplane, which is very similar to the MD-87. The MD-87 airplane is 
powered by two JT8D series turbofan engines, on either side of the fuselage, near 
the tail of the aircraft [28]. A side view of the MD-87 aircraft is shown in Figure 6.1. 

This aircraft would probably be the easiest to place the TCAS antennas on. The 
placement of its engines, between the top and bottom of the fuselage, is in shadow 
of both the top and bottom mounted antennas. Hence, good performance of the 
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Figure 6.1: A scaled sketch of the McDonnell Douglas MD-87 aircraft. Fuselage 
length is shown marked in inches. 


TCAS system may be expected for an antenna placed in the forward region, for 
both top and bottom mounted cases. It is noticed that the wings are slightly above 
the bottom of the fuselage towards the rear of the aircraft and are not expected to 
deteriorate the performance very much. 

6.6.2 Boeing 757-200 

The Boeing 757 series of aircraft is similar to the Boeing 737, in that its two engines 
are placed under the wings. The Boeing 757-200 is a short /medium range transport 
aircraft, 1858" long as shown in Figure 6.2. Note that the fuselage of the Boeing 
737-300 is only 1163" long; whereas, the Boeing 747-200 is 2792" long [28]. 

Due to the increased space available on the top of the fuselage, and placement 
of the VHF blade antenna towards the middle of the aircraft, it is expected that 
the location of antenna for the TCAS system would pose fewer problems than the 
Boeing 737. The bottom antenna on the other hand is restricted by the engine inlets 
and probably should be placed far forward for looking in the nose region. 


104 



Bibliography 


[1] M. T. Pozesky and M. K. Mann, “The US Air Traffic Control System Archi- 
tecture,” Proc. IEEE , vol. 77, pp. 1605-1617, November 1989. Special Issue on 
Air Traffic Control. 

[2] T. S. Perry, “Improving the Worlds’s Largest Most Advanced System,” IEEE 

Spectrum , vol. 28, pp. 22-36, February 1991. Special Report: Air Traffic Con- 
trol. 

[3] R. K. Jurgen, “Antenna Ups Bearing Accuracy in TCASII,” IEEE Spectrum, 
vol. 27, p. 11, April 1990. News Item in ‘Innovations’ Column. 

[4] M. I. Skolnik, Radar Handbook. New York: McGraw-Hill Book Company Inc 

1970. ’’ 

[5] R. G. Rojas, Y. C. Chen, and W. D. Burnside, “Improved Computer Simulation 
of the TCAS III Circular Array Mounted on an Aircraft,” Technical Report 
716199-12, The Ohio State University ElectroScience Laboratory, April 1989. 

[6] J. E. Reed, Enhanced TCAS II Antenna Characteristics Report,” Report 

BCD-TR-090, Allied Bendix Aerospace — Bendix Communications Division 
October 1983. 

[7] J. E. Reed, “Enhanced TCAS II System Summary (A 3D Traffic Alert & Col- 
lision Avoidance System),” Report BCD-TR-098, Allied Bendix Aerospace — 
Bendix Communications Division, April 1984. 

[8] A. I. Sinsky, J. E. Reed, and J. Fee, “Enhanced TCAS II Tracking Accuracy,” 

in AIAA/IEEE Digital Avionics Systems Conference, December 1984 Paoer 
84-2738-cp. V 

[9] K. J. Kunachowicz, “Model Testing of Airborne VHF Direction Finding An- 
tenna System,” in Antennas for Aircraft and Spacecraft, pp. 223-227, IEE Con- 
ference Publication, no. 128, 1975. 

[10] R. G. Kouyoumjian and P. H. Pathak, “A Uniform Geometrical Theory of 
Diffraction for an Edge in a Perfectly Conducting Surface,” Proc. IEEE, vol. 62 
pp. 1448-1461, November 1974. 


106 



[11] J. J. Kim and W. D. Burnside, “Simulation and Analysis of Antennas Radiating 
in a Complex Environment,” IEEE Trans. Antennas and Propagat., vol. AP-34, 
pp. 554-562, April 1986. 

[12] W. D. Burnside, J. J. Kim, B. Grandchamp, R. G. Rojas, and P. Law, “Airborne 
Antenna Radiation Pattern Code User’s Manual,” Technical Report 716199-4, 
The Ohio State University ElectroScience Laboratory, September 1985. 

[13] K. S. Sampath, R. G. Rojas, and W. D. Burnside, “Modelling and Performance 
Analysis of Four and Eight Element TCAS,” Technical Report 716199-15, The 
Ohio State University ElectroScience Laboratory, April 1990. 

[14] R. G. Rojas, K. S. S ampath, and W. D. Burnside, “Analysis and Simulation of 
Collision Avoidance TCAS Antennas mounted on Aircraft,” Applied Computa- 
tional Electromagnetics Society Journal , vol. 6, pp. 112-139, Summer 1991. 

[15] R. G. Rojas and K. S. Sampath, “Analysis of Airframe Effects on the Per- 
formance of TCAS Antennas,” Technical Report 724737-1, The Ohio State 
University ElectroScience Laboratory, August 1992. 

[16] R. F. Harrington, Time Harmonic Electromagnetic Fields . New York: McGraw- 
Hill Electrical Engineering Series, 1961. 

[17] E. H. Newman, “A User’s Manual for the Electromagnetic Surface Patch Code: 
ESP, Version IV,” Technical Report 716199-11, The Ohio State University Elec- 
troScience Laboratory, August 1988. 

[18] A. Altintas, Electromagnetic Scattering from a Class of Open-Ended Waveguide 
Discontinuities. Ph.D. dissertation, The Ohio State University, Dept. Elec. 
Eng., Columbus, OH, March 1986. 

[19] P. H. Pathak, A. Altintas, C. W. Chuang, and S. Barkeshli, “Near Field Scat- 
tering by Rectangular and Circular Inlet Configurations with an Impedance 
Surface Termination,” Technical Report 715267-1, The Ohio State University 
ElectroScience Laboratory, July 1984. 

[20] P. H. Pathak, C. W. Chuang, and M. C. Liang, “Inlet Modeling Studies,” Tech- 
nical Report 717674-1, The Ohio State University ElectroScience Laboratory, 
October 1986. 

[21] C. A. Balanis, Antenna Theory: Analysis and Design . New York: Harper & 
Row Electrical Engineering Series, 1982. 

[22] S. Silver, Microwave Antenna Theory and Design , vol. 12 of Radiation Labora- 
tory Series . New York: McGraw Hill Book Company, Inc., 1949. 


107 



[23] E. K. Sandeman, “Spiral Phase Fields— Possible Applications,” Wireless Engi- 
neer, pp. 96-105, March 1949. 

[24] E. W. Braeden, F. T. Klevenow, E. H. Newman, R. G. Rojas, K. S. Sampath, 
J. T. Scheik, and H. T. Shamansky, “Operation of the Helicopter Antenna Ra- 
diation Prediction Code,” Technical Report, The Ohio State University Elec- 
troScience Laboratory, 1993. To be published. 

[25] P. H. Pathak and N. Wang, “Ray Analysis of Mutual Coupling Between Anten- 
nas on a Convex Surface,” IEEE Trans. Antennas and Propagat., vol. AP-29, 
pp. 911-922, November 1981. 

[26] D. W. Burgess. MIT Lincoln Laboratory, 1991. Personal Communication. 

[27] L. W. Johnson and R. D. Riess, Introduction to Linear Algebra. Reading, MA: 
Addison- Wesley Series in Mathematics, 1981. 

[28] J. W. R. Tailor, ed., Jane’s All the World’s Aircraft. New York: Jane’s Pub- 
lishing Inc., 1988. 

[29] M. E. Peters and E. H. Newman, “A User’s Manual for the Method of Mo- 
ments Aircraft Modeling Code,” Technical Report 716199-14, The Ohio State 
University ElectroScience Laboratory, September 1989. 

[30] J. Van Bladel, Electromagnetic Fields. New York: McGraw-Hill Electrical En- 
gineering Series, 1966. 

[31] R. E. Collin, Foundations for Microwave Engineering. New York: McGraw-Hill 
Electrical Engineering Series, 1966. 

[32] P. H. Pathak and R. J. Burkholder, “Modal, Ray and Beam Techniques for An- 
alyzing the EM Scattering by Open-Ended Waveguide Cavities,” IEEE Trans. 
Antennas and Propagat ., vol. AP-37, pp. 636—647, May 1989. 


108 



Appendix A 

Statistics of TCAS Antennas 


In this appendix, some relevant statistics for the various aircraft are provided. These 
tables provide the location code used in the simulation for future reference. This is 
usually the distance of the antenna from the nose in the UTD model used for the 
aircraft code. The second column is the distance of the TCAS antenna on the actual 
aircraft. Third column gives the angle the antenna normal makes with the forward 
horizontal direction, at the location of the TCAS antenna. The next three columns 
are the output of processing all the data from the aircraft code and they give the 
maximum, average, and the standard deviation of the absolute value of the error for 

each location. 
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Table A.l: Statistics for top mounted CRA TCAS: Boeing 737-300. 



Table A.2: Statistics for top mounted SPA TCAS: Boeing 737-300. 


Location 

Code 


0180 

0210 

0240 

0270 

0300 

0330 

0360 

0390 

0420 

0450 

0480 


Distance 
From Nose 

169.0000 

199.0000 

229.0000 

259.0000 

289.0000 

319.0000 

349.0000 

379.0000 

409.0000 

439.0000 

469.0000 

499.0000 

529.0000 


Theta 

Normal 

99.2469 

97.1508 

95.3603 

93.7608 

92.2791 

90.8633 

89.9599 

89.8543 

89.7483 

89.6418 

89.5346 

89.4263 

89.3167 


Maximum I Absolute I Standard 
[Error | Average Deviation 









Table A.3 j Statistics for bottom mounted CRA TCAS: Boeing 737-300. The model 
for aft mounted antenna used for locations greater than 645" from nose. 


-- 

Location 

Distance 

Theta 

Maximum 

Absolute 

Standard 


Code 

From Nose 

Normal 

| Error | 

Average 

Deviation 


mm 

95.0000 

102.1612 

20.6132 

2.6537 

1.6585 


m 

155.0000 

95.9989 

18.8678 

2.2422 

1.9440 

1 


185.0000 

93.6377 

37.4938 

2.0531 

2.4643 


■B 

215.0000 

91.4692 

17.5247 


2.0902 


mm 

245.0000 

89.9799 

17.0278 




IB 

275.0000 

89.9127 

24.4318 

2.8565 

1 

| 

■si 

305.0000 

89.8453 

22.3935 


IBjRtrfl 


0286 

335.0000 

89.7777 

28.9407 

3.3745 

4.1844 


0316 

365.0000 

89.7098 

24.3487 


mmm 


0346 

395.0000 

89.6414 

32.2942 

4.3415 

6.2711 


0376 

425.0000 

89.5724 

30.1544 

3.6141 



0406 

455.0000 

89.5027 

29.3484 

3.5755 


— 

0441 

490.0000 

89.4202 

24.9285 

3.8472 



0466 

515.0000 

89.3606 

31.9606 

3.9398 

5.9372 


0556 

605.0000 

89.1386 

41.5618 

4.4798 

6.3547 

— 

0321 

645.0000 

88.6294 

27.5043 

4.5429 

6.2296 


0381 

705.0000 

87.8999 

30.8443 

3.7551 

5.9535 


0441 

765.0000 

87.1022 

20.6750 

2.9689 

4.2512 



825.0000 

86.1958 

23.1261 

2.6340 

3.6603 


0561 

885.0000 

85.1138 

20.5303 

2.3502 

2.6607 


0621 

945.0000 

83.7318 

31.8529 

2.5823 

4.4394 


Ill 






Table A.4: Statistics for bottom mounted SPA TCAS: Boeing 737-300. The model 
for aft mounted antenna used for locations greater than 645" from nose. 
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Table A. 5: Statistics for top mounted CRA TCAS: Boeing 747-200. Stations 1380 
and 1440 lie aft of the cupola. 
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Table A. 6: Statistics for top mounted SPA TCAS: Boeing 747-200. Stations 1380 
and 1440 lie aft of the cupola. 


Location 

Station 

Theta 

Maximum 

Absolute 

Standard 

Code 

Number 

Normal 

| Error | 

Average 

Deviation 


270.0000 

115.5353 

5.4403 

4.0316 

1.2404 

ft 

300.0000 

107.4228 

3.7187 

2.4174 

0.8712 


330.0000 

101.8281 

2.8663 

1.5443 

0.6506 


360.0000 

96.1295 

2.1297 

0.8624 

0.4671 

0390 

390.0000 

91.7944 

1.9998 

0.6201 

0.4309 

0420 

420.0000 

93.0957 

2.1131 

0.6440 

0.4604 

0450 

450.0000 

92.2849 

1.9889 

0.6328 

0.4547 

0480 

480.0000 

90.5367 

1.7400 

0.5881 

0.4214 

0510 

510.0000 

89.6165 

1.8281 

0.7403 

0.3830 

0540 

540.0000 

89.3417 

1.8827 

0.6224 

0.3729 

0570 

570.0000 

86.7256 

2.1105 

0.8498 

0.3984 

0600 

600.0000 

85.7099 

2.0913 

0.9906 

0.4379 

0630 

630.0000 

84.2348 

2.5281 

1.0946 

0.3755 

0660 

660.0000 

83.6920 

2.4717 

1.2531 

0.4140 

0690 

690.0000 

82.6091 

2.3826 

1.4119 

0.3658 

0750 

750.0000 

80.5222 

2.9295 

1.7092 

0.3737 

0810 

810.0000 

78.7237 

3.0454 

1.9732 

0.3745 


1380.0000 

90.5189 

32.8940 

3.3374 

5.3732 

fj 

1440.0000 

90.5575 

37.0278 

3.5694 

5.3195 
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Table A. 7: Statistics for bottom mounted CRA TCAS. Boeing 747 


Location 
Code 
0170 
0200 
0230 
0260 
0290 
0320 
0350 
0380 
0410 
0440 
| 0470 

0500 
0530 
0560 
0590 
0620 
0650 
0680 
0710 



Station 

Number 

170.0000 

200.0000 

230.0000 

260.0000 

290.0000 

320.0000 

350.0000 

380.0000 

410.0000 

440.0000 

470.0000 

500.0000 

530.0000 

560.0000 

590.0000 

620.0000 

650.0000 

680.0000 

710.0000 

800.0000 

890.0000 

980.0000 

1070.0000 

1160.0000 


Theta 

Normal 

112.4038 

108.0461 

103.7831 

100.6243 

99.0747 

98.0066 

97.1260 

96.4559 

95.9381 

95.5523 

95.2549 

95.0356 

94.8662 

94.6189 

94.5140 

94.2928 

94.1702 

93.9038 

93.6169 

92.5765 

91.6126 

90.8671 

90.4037 

89.9385 


Maximum 
| Error | 
31.0171 
22.0036 
19.5476 
17.8396 
18.7528 
11.7499 
! 18.3055 

18.4471 
20.2134 
23.1327 
14.2177 
19.5529 
23.0913 
24.3195 
27.0265 
17.1046 
24.2844 
24.6918 
25.5776 
29.4408 


Absolute 

Average 

4.8203 

3.8405 

3.1210 

2.3556 

2.5253 

2.0940 

1.9955 

1.9474 

2.2221 

2.1376 

1.9531 

2.1708 

2.0412 

2.0245 

2.2386 

1.8590 

2.0502 

2.4173 

2.3023 

2.6359 


Standard 
Deviation 
3.1640 
2.3237 
2.4771 
I 1.8371 
2.6292 
1.7301 
2.2995 
2.1145 
2.9699 
2.8784 
2.1094 
2.9524 
3.0873 
3.4849 
3.8630 
2.0605 
3.3240 
4.3770 
4.1770 
5.2018 
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Table A. 8: Statistics for bottom mounted SPA TCAS: Boeing 747-200. 


Location 

Code 

0170 

0200 

0230 

0260 

0290 

0320 

0350 

0380 

0410 

0440 

0470 

0500 

0530 

0560 

0590 

0620 

0650 

0680 

0710 

0800 

0890 

0980 

1070 

1160 


Station 

Number 

170.0000 

200.0000 

230.0000 

260.0000 

290.0000 

320.0000 

350.0000 

380.0000 

410.0000 

440.0000 

470.0000 

500.0000 

530.0000 

560.0000 

590.0000 

620.0000 

650.0000 

680.0000 

710.0000 

800.0000 

890.0000 

980.0000 

1070.0000 

1160.0000 


Theta 

Normal 

”112.4038 

108.0461 

103.7831 

100.6243 

99.0747 

98.0066 

97.1260 

96.4559 

95.9381 

95.5523 

95.2549 

95.0356 

94.8662 

94.6189 

94.5140 

94.2928 

94.1702 

93.9038 

93.6169 

92.5765 

91.6126 

90.8671 

90.4037 

89.9385 


Maximum 
| Error | 
28.8542 
38.9429 
37.2420 
30.1497 
26.2923 
31.9184 
21.6287 
37.3500 
29.5302 
32.1880 
32.2645 
35.6382 
46.7629 
42.7559 
48.3350 
45.4221 
50.0999 
52.8423 
47.5850 
44.9433 
50.5754 
44.6187 
35.8714 
45.9116 


Absolute 
Average 
" 4.7222 
4.1454 
3.5994 
3.1487 
3.0423 
2.6144 
2.6256 
2.9808 
2.9839 
2.8229 
2.9104 
2.8998 
2.8092 
3.2123 
3.9122 
3.4294 
3.3702 
3.2489 
3.3366 
3.8839 
4.3110 
4.4109 
4.6880 
4.7909 


Standard 

Deviation 

3.4279 

3.6746 

3.9984 

3.9460 

3.5636 

3.1870 

3.5423 

4.7845 

4.3333 

4.3041 

4.3667 

4.0981 

4.8430 

6.0199 

7.6843 

6.5876 

6.4994 

6.2358 

6.1165 

7.5809 

7.5467 

7.1810 

7.4140 

8.0948 
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Table A. 9: Statistics for top 


mounted CRA TCAS: Boeing 727-200. 


Location 

Distance 

Theta 

Code 

From Nose 

Normal 

0100 

100.0000 

102.3022] 

0140 

140.0000 

98.9772 

0180 

180.0000 

96.5303 

0200 

200.0000 

95.4916 

0220 

220.0000 

94.5344 

0260 

260.0000 

92.7874 

0300 

300.0000 

91.1730 

0320 

320.0000 

90.3896 

0340 

340.0000 

89.9625 

0380 

380.0000 

89.8125 

0420 

420.0000 

89.6616 

0460 

460.0000 

89.5094 

0500 

500.0000 

89.3549 

0540 

540.0000 

89.1976 

0580 

580.0000 

89.0366 


Maximum 

|Error| 

Absolute 

Average 

Standard 

Deviation 

3.9868 

1.5764 

0.5341 

2.4316 

1.1656 

0.4360 

2.4374 

0.8667 

0.5235 

2.3658 

0.9394 

0.4318 

2.1595 

0.8669 

0.4356 

2.9008 

0.9677 

0.4569 

3.6957 

1.0707 

0.6966 

3.9930 

1.0916 

0.6490 

3.7680 

1.2352 

0.7462 

4.0563 

1.4472 

0.7959 

4.7993 

1.4166 

0.8804 

4.9945 

1.4267 

0.9297 

5.5235 

1.5664 

1.0979 

6.5893 

1.6835 

1.2862 

7.0991 

1.9211 

1.4842 


Table A. 10: Statistic, for top mounted SPA TCAS: Boeing 727-200. 


Location 

Code 

Distance 
From Nose 

Theta 

Normal 

Maximum 
| Error | 

Absolute 

Average 

Standard 
Deviation 
i oocn 

0100 

0140 

0180 

0200 

0220 

0260 

0300 

0320 

0340 

0380 

0420 

0460 

0500 

0540 

0580 

100.0000 

140.0000 

180.0000 
200.0000 
220.0000 
260.0000 

300.0000 

320.0000 

340.0000 

380.0000 

420.0000 

460.0000 

500.0000 

540.0000 

580.0000 

102.3002 

98.9761 

96.5298 

95.4911 

94.5340 

92.7872 

91.1729 

90.3896 

89.9625 

89.8125 

89.6617 

89.5094 

89.3549 

89.1976 

89.0367 

15.0834 

6.7901 

8.3212 

8.9065 

12.1936 

14.0352 

15.4140 

20.7349 

17.5436 

20.8014 

24.6850 

22.0273 

23.6262 

32.5785 

41.5553 

2.1480 

1:8592 

1.9102 

2.0151 

2.2136 

2.5460 

3.0154 

3.6282 

4.3637 

4.0296 

4.6279 

4.7512 

4.8537 

6.1714 

7.3432 

1.0037 

1.3128 

1.7084 

2.2104 

2.8560 

3.5035 

4.5359 

4.2194 

4.8206 

6.0154 

5.5725 

5.6560 

7.6192 

9.4367 
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Table A. 11: Statistics for bottom mounted CRA TCAS: Boeing 727-200. 


Location 

Code 

Distance 
From Nose 

Theta 

Normal 

Maximum 
| Error | 

Absolute 

Average 

Standard 

Deviation 

0100 

100.0000 

100.0279 

45.0644 

2.9680 

5.6467 

0140 

140.0000 

97.2995 

40.0294 

2.7964 

6.2837 

0180 

180.0000 

95.3030 

58.3248 

3.0624 

7.9488 

0220 

220.0000 

93.6793 

24.8572 

1.2885 

1.7567 

0240 

240.0000 

92.9514 

8.0948 

1.1278 

1.0205 

0260 

260.0000 

92.2608 

14.8468 

1.0957 

1.1760 

0300 

' 300.0000 

90.9512 

9.0197 

1.0902 

1.1176 

0340 

340.0000 

89.9696 

24.4167 

1.2039 

1.9215 

0360 

360.0000 

89.9088 

12.6119 

1.1121 

1.2853 

0380 

380.0000 

89.8480 

30.9790 

1.6472 

3.4001 

0420 

420.0000 

89.7257 

33.1766 

2.7763 

5.7794 

0460 

460.0000 

89.6022 

69.0273 

3.2817 

8.0215 

0480 

480.0000 

89.5398 

90.8555 

7.4584 

14.3458 


Table A. 12: Statistics for bottom mounted SPA TCAS: Boeing 727-200. 


Location 
Code 
" 0100 
0140 
0180 
0220 
0240 
0260 
0300 
0340 
0360 
0380 
0420 
0460 
0480 


Distance 
From Nose 
100.0000 

140.0000 

180.0000 
220.0000 

240.0000 

260.0000 

300.0000 

340.0000 

360.0000 

380.0000 

420.0000 

460.0000 

480.0000 


Theta 

Normal 

"100.0260 

97.2985 

95.3024 

93.6789 

92.9511 
92.2605 

90.9511 
89.9696 
89.9088 
89.8480 
89.7257 
89.6022 
89.5398 


Maximum 
| Error | 
34.9310 
24.9701 
39.3768 
31.0480 
26.1150 
33.9125 
31.9772 
37.6930 
32.4067 
28.1335 
55.0527 
68.5810 
70.0208 


Absolute 

Average 

2.8616 

2.6363 

2.7321 

1.8654 

1.3524 

1.8206 

1.6031 

2.4063 

1.4347 

2.0772 

3.1774 

3.8541 

6.6051 


Standard 
Deviation 
" 3.8932 
3.9149 
4.8413 
2.3767 
2.2934 
2.5723 
1.7913 
3.0135 
3.0853 
2.7737 
6.1373 
7.1191 
12.6037 I 
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Appendix B 


Helicopter Antenna Radiation 
Prediction (HARP) 

B.l Introduction 

The Helicopter Antenna Radiation Prediction (HARP) is the name given to a com- 
puter code developed at the Ohio State University ElectroScience Laboratory [24]. 
It integrates the techniques of modern computer graphics with two electromagnetic 
far zone scattering and radiation analysis programs. The result is an easy to use, 
versatile, package to graphically complement and help visualize the mathematical 
models and patterns at various stages of the program. 

Two widely used computational techniques of electromagnetic analysis are the 
method of moments (MM) and the uniform geometrical theory of diffraction (UTD). 
MM is a low frequency method, used when the objects are not electrically large. 
For the sources radiating in the vicinity of electrically large bodies, UTD is more 
appropriate. For a given problem, hence, the method chosen depends upon the size 
of the body as well as available computational resources. 

HARP is an X windows based graphical user interface (GUI) that enables a user 
to build a model of an aircraft from the blue prints of the cross sections of the body 
along its length. The user can then predict the electromagnetic radiation patterns 
of specified antennas on the airplane or helicopter. Though HARP was developed 


119 



THE HELICOPTER ANTENNA RADIATION PREDICTION SYSTEM 



Figure B.l: Block diagram of the major components of HARP. 


primarily for modeling helicopters, the modeling methods are very general and hence 
applicable to many varieties of bodies, of essentially arbitrary shape, like airplanes, 
missiles etc. 


B.2 Overview 

Figure B.l depicts a block diagram of HARP. The major functions are denoted by 
their own blocks. Most of the modeling involves the entry and editing of the body 
data until a satisfactory body is built. Note that the only electromagnetics codes 
are blocks, labeled Perform MM Computations and Perform UTD Computations. 
There are three main parts in the block diagram given in Figure B.l: 
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1. User enters the aircraft geometry as cross sections of the body along some 
sections. Fins, wings, stabilizers, etc., are entered as flat plates. A continuous 
surface representation of the fuselage is generated by using splined interpola- 
tion of the cross sections. This surface is then displayed on the workstation 
screen along with the flat plates. HARP does not yet have the capability to 
read the various databases used by the aviation/aerospace industry, but is 
anticipated soon. 

2. Depending on the frequency and size, one uses either the high frequency 
method (UTD) or the low frequency MM technique. The models are different 
based on the method chosen, as given below: 

(a) If the aircraft is to be analyzed via UTD, a model that simulates the 
fuselage by a composite ellipsoid is generated and displayed. The radia- 
tion patterns of the antennas on this model are then computed using the 
airborne antenna code [12]. 

(b) If the user chooses MM, then the aircraft is modeled by polygonal plates 
and displayed. The radiation patterns are then computed via MM using 
the Aircraft Modeling Code (AMC) [29]. 

3. The results of the EM computation are displayed in a user friendly form, that 
is both easy to interpret and analyze. 

B.3 Low Frequency MM Analysis 

When conditions dictate that MM be the method most suitable for the analysis, 
the first step is the construction of a model to be used with AMC. This program 
requires that the fuselage be defined by cross sections, wings and fins by flat plates 
and antennas, rotor, etc. by wires. HARP converts the continuous surface model 
into a set of cross sections compatible with AMC. This is not a trivial task because 
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AMO needs the fin and wing attachments to exactly match the grid formed by the 
fuselage. AMC then converts the fuselage cross section data into a set of polygonal 
flat plates. The wings and fins are also flat plates. The rotor is modeled by a set of 
arbitrary wires. 

Once the surface patch model of the helicopter has been computed, AMC con- 
tinues with the complete MM computation. This involves the determination of the 
impedance matrix, the right hand side voltage vector, as well as the matrix of solu- 
tion for the patch and wire currents. Once this is computed, AMC finally computes 
the required patterns in the far zone. Depending on the electrical size of the body, 
and the computational power at hand, the MM computation can take a few seconds 
to several hours of execution time. 

B.4 High Frequency UTD Analysis 

If the body is electrically large, then UTD provides the solution more efficiently. 
This task is accomplished by the use of NEWAIR3 computer program. NEWAIR3 
requires the fuselage to be modeled by a composite ellipsoid and the remainder of 
the structures by flat plates. The composite ellipsoid need accurately model the 
aircraft only the neighborhood of the antenna. Note that the antenna is required 
to be attached to the fuselage. Flat plates (if necessary) need only be approximate. 
Both attached as well as unattached plates to the fuselage are treated by NEWAIR3. 

It is emphasized that, in contrast to the moment method technique where the 
entire model has to be reasonably accurate, the UTD model needs to be accurate 
only in the vicinity of the antenna. This is due to the more localized nature of high 
frequency phenomena, such as diffraction, reflection etc. Therefore, the airplane 
model can be less accurate in the shadow regions, and locations away from the 
antenna. 

HARP generates this composite ellipsoid very accurately, as well as the flat plates 
that are necessary to reasonably model the aircraft. If the user is satisfied with 
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the views of the model presented on the screen, he can then continue on to the 
computation of the required radiation patterns. These can once again be displayed 

on the screen. 
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Appendix C 

Scattering from Waveguide Inlets 


The analysis of the near field bistatic EM scattering by the waveguide is primarily 
based on the GTD and the ECM, as well as the PTD, combined with the self- 
consistent MSM. In the MSM based analysis, the basic scattering mechanisms are 
isolated and identified as being associated with the scattering by junctions I and II 
as in the waveguide shown in Figure C.l. The MSM procedure requires knowledge 
of the generalized scattering matrices [5„], [S ]2 ], [S 21 ], [S 22 ], and [S, ] for the two 
junctions. 

The scattering matrix concept in microwave circuit analysis is based on the prop- 
agating modes within the circuit. It characterizes the circuit properties or waveguide 
discontinuities of a microwave network. Hence, the circuit scattering matrix for mi- 
crowave networks is limited to guided wave (interior) regions. On the other hand, the 
polarization scattering matrix is defined for exterior regions [30]. The generalized 
scattering matrix extends the concept of the microwave circuit scattering matrix 
to include evanescent modes for interior regions as well as the polarization scat- 
tering matrix for external scattering. Thus, this approach is useful for problems 
that involve a coupling between interior and exterior regions as in the geometry of 
Figure 5.1. 
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Figure C.l: Illustration of multiple scattering between junctions I and II in a term)- 
nated waveguide. 

Let E' denote the electric field intensity of the incident EM plane wave. One 



may express E 1 as: 

E* = 6E' 0 + <£e; 


(C.l) 

— 

with 

E' 0 = A e e- }kif , 

and 

(0.2) 

— 


e; = A^ k ' f 


(C.3) 

— 

where 

V 

= — fc(sin 9' cos fix + sin 9 * sin 

<j> % y -f cos 8'z) 

(0.4) 

— 

and 

f = xx + yy + zz. 


(C.5) 


Note that k is the free space wave number and f is the position vector of an obser- 
vat ion point at P. 

The scattering matrix [S u ] is described next. The matrix [Sji] relates the electric 
field E*° scattered from only the open end at 2 = 0 to the field E* which is incident 
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Figure C.2: Scattering by only the rim at the open end of the inlet. This scattering 
affects both the exterior and interior regions. 

upon the open end as follows (see Figure C.2): 


' E? ' 


' S% 

CSO * 
^x* 


r An 

K 


C so 
D y0 

Cso 

°v4> 


Ai 

. E z° . 


O so 

L ° z$ 

Cso 

^ z 4> J 




where 

E so (P) = E 9 x °x + Ey°y + E 9 z °z 


(C.6) 


(C.7) 


and 



' G,«o 
°x0 

Cso " 

°x<j> 

[‘S’ll] - 

Cso 

°y0 

Cso 


1 — 

Co 

Qao 

Oao 

^ * 4 > J 


(C.8) 


Note that Ag and A $ are the amplitudes of the incident fields Eg and Ej,, respectively. 

The scattering matrix [5* j 2 ] converts the waveguide modal fields incident at the 
open end (2 = 0) from the interior region ( 2 < 0) to the fields radiated by these 
modes from the open end as shown in Figure C.3. The modal electric field E* within 
the waveguide region (2 < 0) may be represented as [31]: 


E g = HH\ A tm^ n m^ i0 ' nmZ + BL(e»m ± ze* nm )e* iflnmt ] 


(C.9) 
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Figure C.3: Radiation and reflection at the open end due to a waveguide mo 
incident on the opening. 

where, e'„ m end e„„ denote the transverse (to z) electric vector mode functions for 
the TEnm and TM„„, modes, respectively. The Ai„ and in (C.9) are the mo a 
amplitude coefficients, and ft. and ft. denote the propagation constants for the 
TE nm and TM„„ modes, respectively. Here, » and m refer to the modal indices 

associated with the transverse eigenvalues in « and , and eigenfunctions (modes), 

. . . , nf<k TM modal electric field. The corresponding 

and e zn m is the z-component ot the lM„ m mou. 

e' for the TE nm are identically zero by definition. The superscripts + and - 
above refer to the modes propagating in the +S and -i directions, respectively. 

It E'"(P) denotes the electric field at P external to the waveguide region, which 
is radiated by the modes striking the open end, then the Scattering matrix |S„] 

relates E r "(P) to EJ by 

r is;"'] [s;l ’ 

= is;"'] |s;l 
is;"'] is,™] J 


b;° 

K 

E’° 


l wj 1 


(C.10) 


where 


E ro (P) = E r x °x + E[°y + E r z °z 


(C.ll) 
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and 


with 


and 



Mul = 

Mrl 

Mr) 1 



Mr'i 

m;i 




Mri 

Mr] J 


[sr’i 

= f^iOl $102 ' 

•• s;„ 

sin •• 

■ s;„„ ••■] 

mti 

= [‘S'tOl 5 'i02 • 

* * Sin 

5"ii2 • ' • 

Sinm * ’ •] 


(C.12) 


(C.13) 


(C.14) 


where i = x, y, z. 

The scattering matrix Mr,] describes the transformation or the coupling of the 
incident plane wave field into the waveguide modes as shown in Figure C.2. It is 
also noted that the determination of |S„] is related to the determination of [S,,] via 
reciprocity, when the receiver recedes to the far zone in the latter problem. It may 
be seen that [S 21 ] relates E~ to E' via 


' M;J 1 


M„-„] . 



[SL.] 

l^enm] 


[5, 

& 


*^nm] I | Ag 

'<£nm] J [ A,], 


(0-15) 


where 


IS,.] = f Mft,„] M*,J 
' 1 MU MU 


(C.16) 

The scattering matrix [S 22 ] is a modal reflection coefficient matrix associated with 
the mechanisms depicted in Figure C.3. In particular, the elements of [S 22 ] describe 
the reflection coefficients associated with the TE„ m and TM nm modes reflected back 
from the open end (at z = 0) into the waveguide region (z < 0) when either a TE W 

or a TM m mode in incident on the open end from within the waveguide. Thus, the 
matrix [5 22 ] relates E~ to E + bv 

3 g J 
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where 


[£22] = 


MK.J MU 
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*hr 1 
nm;pq\ 

nm;pq\ 


(C.17) 


(C.18) 
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Figure C.4: Reflection of waveguide inodes from the termination at 2 L 


Note that the reflection coefficient R h n c m;pq denotes the following: A TM W (or e 

type) mode with modal amplitude B+ which is incident at the open end is partly 
transformed (coupled) into a reflected TE nm (or h type) mode with modal amplitude 


~''nm;pq pq 

The scattering matrix [Sr], like [S»J, is also a reflection coefficient type matrix 
which is associated with the termination at z = -L. This mechanism is lUustrated 


in Figure C.4. For this case, one obtains that 



Wl 

1^] J 


(C.19) 


where 


[Sr] 


' [r“ ,„] [r^J' 


As before, for example, 


r 


he 

nm\pq 


lates A+ to B m via 


(C.20) 


At any given frequency, the waveguide region (interior) can support a finite num- 
ber of propagating modes and an infinite number of non-propagating (evanescent) 
modes. Hence, strictly speaking, the matrices (S 12 ], l&i], [^ 22 ] and [S\ ) are of infi 
nite order. However, in practice, one need retain only a finite number of elements 
of the matrix (finite number of modes), because the distance L is generally large, 
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E mo — — 
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E 21 
~ E 12 


INTERIOR 
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[P]E 21 


JUNCTION I 

z=0 


JUNCTION II 

z=~L 


Figure C.5: Multiple interations between junctions I and II for calculating the total 
scattered field. 


and the infinite number of modes generated at junction I contribute negligibly at 
junction II, and vice versa. Therefore, only a finite number of propagating modes 
are retained in practice. If the distance L is too small, such that these evanescent 
modes do become important, then one must include these modes, but one may still 
ignore the higher-order evanescent mode contributions. 


The fields scattered by the inlet may be viewed as the superposition of two parts, 
one of which is the scattered field from the open end (rim) by itself. The other field is 
that radiated by the open end as a result of all the interactions between the open end 
and the termination, which arises from the incident field coupled into the waveguide 

region. Referring to Figure C.5, it i6 seen that the field scattered by junction I is 
simply 


[E“] = [$„][£•'] 


where 


[E ao ] = 


E*° 

K 

E‘° 


(C.21) 


(C.22) 


A part of the incident field is coupled at I into the waveguide region, which 
becomes the incident field at junction II from which it is reflected. The field reflected 
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by junction II impinges again on T. where it experiences further scattering into the 
exterior and interior regions. This process repeats ad infinitum. The Reids resulting 
from these multiple interactions may be found from the self-consistent MSM analysis, 
which is briefly described below. 

Let [Eu] be the total field incident at I from II after taking all the multiple 
interactions into account. Similarly, let [ft,] represent the total field incident at 
junction II from I after taking all the multiple interactions into account. Then, 
the total scattered field in the exterior region [ft] consists of superposition of [ft»] 
and [ft" 0 ], where [£ mo ] is the field scattered into the exterior region when ( ft, ! 

incident on I. Hence, 

\E 3 ) = [E*°] + [ E mo ] (C.23) 


where 


[E mn ] = [Si 2 ][£i 2 ]| 2=o 
W riting the expressions for [E l2 ] and \E- n } in terms of [S 2 1 ] and [S 22 ], 

\E 2 i]| = [S-n}[E'\ + [5 22 ]l^. 2 ]| 8=o 


(C.24) 


(C.25) 


and 

[£r,x]| _ o = [P][s, ][/’]lft.l |„ 0 (c ' 26) 

where, [Sr] denotes the reflection coefficient scattering matrix for the junction II as 
defined in (C.20), and [P| is a diagonal matrix accounting for the phase delay along 
the length L of the waveguide between junction. I and II. Eliminating [ftrl from 
(C.25) and (C.26) yields: 

(C.27) 


((/] - [SaMScm [ft,] = „ = 


where [/] is an identity matrix of infinite order. Solving (C.27) for [ft,], one obtains 
that 

[ft,l| =([/]-[ftr][P][ft]|P]r' [ft, ][£']■ < c - 28) 


2=0 
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Substituting (0.28) into (0.2fi), one finds that 




z=0 


[^][5r][P]([/] - ^[PJ^rKP])- 1 ^!]^]. 


(C.29) 


From (C.29) and (C.24), the total exterior scattered field from the internal scattering 
is given by 


[£"“] = PuMSrllPKW - [52,][P][Sr][/’])- | [S ! ,][F]. (C.30) 

Finally, the self-consistent expression for the total scattered field [£'| can be 
written from (C.21), (C.30) and (C.23) as 


[ £ *1 = {[Sn] + [S.iMSrMM - [Sn|[P](Sr][/ , ])‘ , (S ; ,]} [£']. (C.31) 

It is noted that the above scattering matrices are different from the conventional 
scattering matrices due to the fact that they contain range information. They remain 
valid for field calculations in the near zone of the inlet opening. Further, as noted 
in [20], the above analysis is restricted to angles of incidence and scattering which 
are about 60 0 or less from the forward axis of the inlets. 

The scattering matrices [5u], [ 512 ], [S^i], and [ 522 ] can be found via asymp- 
totic high frequency techniques (such as UTD, ECM and PTD). The methods are 
discussed in detail in [32] and [20] and not repeated here. 
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